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PREFACE 


The  High-Speed  Research  Program  sponsored  the  NASA  High-Speed  Research  Program 
Aerodynamic  Performance  Review  on  February  9-13,  1998  in  Los  Angeles,  California. 
The  review  was  designed  to  bring  together  NASA  and  industry  High-Speed  Civil 
Transport  (HSCT)  Aerodynamic  Performance  technology  development  participants  in  areas 
of:  Configuration  Aerodynamics  (transonic  and  supersonic  cruise  drag  prediction  and 
minimization),  High-Lift,  and  Flight  Controls.  The  review  objectives  were  to:  (1)  report 
the  progress  and  status  of  HSCT  aerodynamic  performance  technology  development;  (2) 
disseminate  this  technology  within  the  appropriate  technical  communities;  and  (3)  promote 
synergy  among  the  scientist  and  engineers  working  HSCT  aerodynamics.  In  particular, 
single-  and  multi-point  optimized  HSCT  configurations,  HSCT  high-lift  system 
performance  predictions,  and  HSCT  Motion  Simulator  results  were  presented  along  with 
executive  summaries  for  all  the  Aerodynamic  Performance  technology  areas.  The  HSR  AP 
Technical  Review  was  held  simultaneously  with  the  annual  review  of  the  following 
airframe  technology  areas:  Materials  and  Structures,  Environmental  Impact,  FlightDeck, 
and  Technology  Integration.  Thus,  a fourth  objective  of  the  Review  was  to  promote 
synergy  between  the  Aerodynamic  Performance  technology  area  and  the  other  technology 
areas  within  the  airframe  element  of  the  HSR  Program. 

The  workshop  was  organized  in  three  sections  as  follows: 

Section  I Independent  Sessions 

Section  II  Plenary  Session 

Section  HI  Executive  Summaries 

The  work  performed  in  the  Configuration  Aerodynamics  element  of  the  High-Speed 
Research  Program  during  1997  was  presented  in  the  following  sessions: 

Analysis  Methods  and  CFD  Validation 
Viscous  Drag  Predictions  and  Testing  Methods 
Aerodynamic  Design  Optimization  Capability 
Nacelle/Diverter  Design  and  Airplane  Integration 
Configuration  Assessments  and  Fundamental  Studies 
Technology  Integration  (TI)  Studies  related  to  Configuration  Aerodynamics 
(CA  / TI  Joint  Session) 

The  work  performed  in  the  High  Lift  (HL)  element  of  the  High-Speed  Research  Program 
during  1997  was  presented  in  the  following  sessions: 

Concept  Development 
Test  Programs  and  Techniques 
Analytical  Methods 
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The  proceedings  for  the  Aerodynamic  Performance  Annual  Review  are  published  in  two 
volumes: 

Volume  I,  Parts  1 and  2 Configuration  Aerodynamics 

Volume  II  High  Lift 

AP  Review  Chairperson:  Naomi  McMillin 

NASA  Langley  Research  Center 
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Time  accurate  linear  and  Euler  dynamic  ground  effect  analyses 
PD  method  (AER02S  based)  not  reliable  - TCA  not  well  predicted 
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Concept  Development  - Outboard  Sealed  Slat 
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Navier-Stokes  provides  good 
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Aft  fuselage  balance  to  directly  measure  tail  loads 
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Pre-test  analysis  of  TCA2.8-28  configuration  has  been  conducted  to  assess  the  effect  of 
planform  variation  on  the  high  lift  performance  and  to  support  the  upcoming  TCA-4  test. 
The  geometric  difference  between  the  TCA2.8-28  and  the  baseline  TCA  occurs  in  the 
outboard  wing  panel  with  a lower  leading  edge  sweep  and  a higher  aspect  ratio. 
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Effects  of  Flaps  on  TCA  Flow  Field 
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Lower  Surface  Limiting  Streamlines  for  TCA 

Configuration 
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Tap  Locations  for  the  TCA  Configuration 
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Chordwise  Pressure  Comparison  for  TCA  30/10  (Part  Span) 
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Chordwise  Pressure  Comparison  for  TCA  30/1 0 (Part  Span) 

HSCT  Aerodynamics,  Long  Beach 
M=0.3,  AoA=10  deg,  Re=8  Million  (MAC) 


/ r~i 


o 

CNJ 

II 


M 


m o 
ox 


O 

LO 

fl 

> 


in  o 
ox 


Spanwise  Pressure  Comparison  for  TCA  30/10  (Part  Span) 
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Spanwise  Pressure  Comparison  for  TCA  30/1 0 (Part  Span) 

HSCT  Aerodynamics,  Long  Beach 
M=0.3,  AoA=10  deg,  Re=8  Million  (MAC) 
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Effect  of  LE  Flap  Angle  on  Surface  Pressure 
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Spanwise  Extent  of  the  LE  Flaps  for  TCA  Configurations 


c 

0 

0 

-Q 


C 0 

|lS8 

0)  C ^ 

5 a® 

O 3 T- 

£ 05  8 

o n ^ 

III 

0 0)  H- 

(0  0)0 
9-u  c 
« <D  g 

**“■  O)  V 

0 C O 
CD~  © 

-D  2 h- 
<D  ® O 

Z1  CD  *D 
CD  — 

fit 

1 

— tr  £ 

® o Z 

£ Cl  ca 
o 03  c 

0 J2  2 

"5<  CD 

® •o  .g> 

o **= 
$ E § 

c c < 
0 O o 
C1‘5  h 

•H  0 

S 3.1 

5 I £ 

■«  O -o 
© • 52 

!t  35  CO 
CD  ^ C 

£ a © 

I-  m ra 


1956 


o> 

c 

o 

0)  -J 

S CO 

O o 

* 1 

5 I 

O)  -o 


1 


Effects  of  the  LE  Flap  Extent  on  the  Vortical  Flow 
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Effect  of  Inboard  LE  Flaps  on  the  Flow  Characteristics 
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Effect  of  Spanwise  Extent  of  LE  Flaps  on  Force  and  Moment 
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Effect  of  Spanwise  Extent  of  LE  Flaps  on  Drag  and 

Performance 
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Inboard  LE  Droop/Camber  Variation  for  Performance  Study 
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Inboard  LE  Droop/Camber  Variation  for  Performance  Study 

HSCT  Aerodynamics,  Long  Beach 
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TCA  Inboard  LE  Modifications 
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It  is  noted  thet  the  chordwise  extent  of  the  modified  region  for  the  cembered  configuration 
is  significantly  larger  than  that  of  the  dropped  configuration.  Therefore,  this  study  will  also 
provide  some  insight  for  the  future  study  on  the  effect  of  the  chordwise  extent  of  the 
leading  edge  flaps  on  the  flow. 
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Effects  of  Inboard  Droop/Camber  on  Flow  Separation 
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Effects  of  Inboard  Droop/Camber  on  Flow  Separation 

HSCT  Aerodynamics,  Long  Beach 

TCA  30/10  Configurations 
M=0.3,  a=12  deg,  Re=8  million  (MAC) 
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Effects  of  Inboard  Droop  Camber  on  TCA  Flow  Solutions 
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Chordwise  Pressure  Variation  due  to  Inboard  Droop/Camber 
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Chordwise  Pressure  Variation  due  to  Inboard  Droop/Camber 

HSCT  Aerodynamics,  Long  Beach 
M=0.3,  a=10  deg,  Re=8  million 
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TCA2.8-28  Flow  Analysis 
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Planform  Comparison  for  the  TCA  Configurations 
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Planform  Comparison  for  the  TCA  Configurations 

HSCT  Aerodynamics,  Long  Beach 
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AER02S  Flap  Optimization  Results  for  TCA2.8-28 
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Computational  Surface  Grid  for  the  TCA2.8-28  25/1 0 

Configuration 
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Flow  Solution  for  the  TCA2.8-28  25/00  Configuration 
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TCA2.8-28  25/10  Flow  Solutions  at  a-sweep 
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Lift  & Drag  Comparison  for  TCA2.8-28 
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separated,  and  the  discrepancy  magnifies  at  higher  angles  of  attack. 


Lift  & Drag  Comparison  for  TCA2.8-28 
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Add  CD0=93  counts  to  AER02S 


Pitching  Moment  comparison  for  TCA2.8-28 
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Pitching  Moment  Comparison  for  TCA2.8-28 
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TCA2.8-28  Outboard  TE  Flap  Characteristics 


ID 

CD  ^ 

O UJ 

■§  m 

o d. 

1 ^ °- 

§£ 

» 2 ^ 

■Q  o 

w i: 

5 ^ 

c ~ 

o co 

„ 4- 

4— ' _Q 

O o 

<D  ° 
3=  C 
<D  0 
TJ  0 
0.-0 
0 0 

0 ® 
O)  -o 

' -O  C0 
<D  C 

CD  .2 


'5  0 

i| 

CO  4- 

s * 

4— • 

S'" 

a)  c 
5 .2 
■2  -g 

> CT) 


0 

in  £ co 
t—  £ to  co 
-O  o E O) 
c o o .2 

(u  **“"  4— » 4— » 

_ o co  o 
o u a m 

'T  Q.  $ Q- 
co  « o i5 

« O m O 
n CT  n CT 
m XJ  co  TD 
^ O t ® 
CD  CD  £ CD 
CD.£  W .£ 
TJ  = ir  = 
CD  CO  ® « 

O)  •*"’  Q.  ■*"J 

c u c p 


C O 
.2  ° 
co  CM 

|<° 

.E  CM 

CO  < 

*co  P 
o I— 


1 § 

JD 

a £ 

CO  o 

E E 
o r 

O o 
as  c 
£ o 

O)  O 

c c 

2 co 
£ co 

II 

CO  CO 

<D  c 
' 0 0 

cd  _r 

CD  £4 
*D  = 

O CO 

^ O 
LO 

CM  CD 


CD  CO 

-C  o 

i-  a 

. 3 

CD  O 

o 0 

.2>£ 

*+- 

CD  O 


O H 
c 

O CO 
~ CD 
O o 
03  ^ 
3=  o> 

CD  CD 
*D  "O 


.2  3 

c ° 

I - 

O Q) 

« I 

CD  — . 

« c ' 

CO  CO  ■ 
CD  »_ 

0)  CD  : 

o)  §-• 

CD  £- 
T5  3 1 

in  £ 

<1~""  4— » 

■D  O . 
C *-  ' 
CO  O) 

o .2 

i-  T3 

II  C 
CO  O 
_£Z  CL 
CL  CO 
r=  CD 
CO  vC 


0)  <d  co 

•g  .=  ® 
« ~ o> 

i=5  5 c 

CO  ° CO 

o "2  c 

E _ .8  .2 

o 2 « t5 
^ (8  ® 
C £ Q.  5= 

CD  O CD  0> 
0 ° CO  "□ 

W ® (o  a 
to  ;C  © CO 

& § - 

0 0 E .23 
■O  g3  5 ^ 

c m cd  o) 

-—  T—  <—  Q. 

m II  — O 
2 to  (0  M 

*“*  +z  CO  J. 
"E  10  co  IE 
CO  "O  Q) 

0 0 « 2 

1 « §•£ 

^ 8 ® 2 

- CO  g "O 

C CD  2 r- 

d ° e “ <o 

0 £ O 3 C0 

+=  m O (fl  — 

CO  S 0 »_  0 

1 IIP. 

~ M—  4-4 


£ * &3 

m r iS  to 

£ p ® l- 
.era 

& <o  r « 
®g°< 

O)  0 to  CD 
^ CT  03  O 
73  0 £ aj 

O -O  ® *t 

*1 1 co  s 2 
JL  II  o w 
S to  ® a 

+;  +-  t Cd 

<003= 


° 8 


.§  2 

w i: 

< O 

• "D 
03  c 
3 

T3  o 

0 _Q 

05  m 
■2  | 
■a  > 
0 ^ 

^ -o 


0 3 

T5 

0 += 


1992 


deflections  implies  that  the  flap  effectiveness  due  to  the  upper  surface  vanishes  when  the 
outboard  wing  is  completely  separated.  The  additional  flap  control  power  at  a higher 
deflection  angle  is  primarily  attributed  to  the  higher  pressure  acting  on  the  flap  lower 
surface  and  varies  linearly  with  the  deflection  angle. 
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HSCT  Aerodynamics,  Long  Beach 

M=0.3,  Re=8  million,  8(LE/TE)=25/1 0 deg 
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O/B  TE  Flap  Deflection  Angle  (degrees) 


Baseline  TCA  Outboard  TE  Flap  Characteristics 
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Baseline  TCA  Outboard  TE  Flap  Characteristics 

HSCT  Aerodynamics,  Long  Beach 
M=0.3,  Re=8  million,  8(LE/TE)=30/1 0 deg 
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Surface  Flow  Patterns  at  an  Off-Design  Condition 
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Surface  Flow  Patterns  at  an  Off-Design  Condition 

„ „ *—  ■ HSCT  Aerodynamics,  Long  Beach 

M=0.3,  AoA=15  deg,  Re=8  million 
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f(degree  of  camber,  chordwise  extent,  curvature,  LE  radius) 
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V-tail  Area  415  sq.  ft. 

Canard  220  sq.  ft. 

Chin  fin  25  sq.  ft. 
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In  the  case  of  the  PTC  buildup,  the  Ref  H baseline  data  set  only  includes  the  LE/TE  flap  combinations  of  00/00, 
00/30,  and  30/10.  Any  other  flap  settings  for  the  PTC  are  predicted  using  AER02s  to  increment  from  one  of 
these  available  baseline  settings. 

In  the  case  of  the  programmed  flap,  the  process  is  used  169  times  (13  leading  edges  and  13  trailing  edges)  to 
create  169  different  polars  that  are  interpolated  to  find  the  LE/TE  combination  giving  the  best  L/D. 
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The  dot-dash  line  represents  the  trimmed  polar  if  there  were  no  canard.  This  data,  however,  is  somewhat 
pessimistic  because  the  tail  size  on  the  PTC  would  not  be  large  enough  without  the  canard.  A larger  tail  would 
improve  the  trim  drag  characteristics  of  the  tail-only  airplane. 
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169  polars  are  created  separately  by  running  through  the  buildup  process  mentioned  previously.  Each  of  these 
polars  is  then  interpolated  at  a series  of  lift-coefficients  to  determine  which  flap  combination  offers  the  best  L/D 
at  each  lift-coefficient.  The  resulting  polar  is  termed  the  “status”  programmed  flap  polar. 

At  each  of  the  lift-coefficients  the  angle-of-attack  for  the  best  flap  combination  is  also  determined. 
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Correlation  Of  CFD  Calculations 
And  Wind  Tunnel  Measurements 
For  The  M2.4-7A  Arrow  Wing  Configuration 
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HIGH  LIFT  TECHNOLOGY  DEVELOPMENT  (TASK  33) 

This  study  falls  under  the  subject  areas: 


CO 

~o 

o 


a> 

E 


CO 

CD 


CO 

CD 

“> 

03 


CO 


2057 


2058 


A0ST/AV 


OUTLINE 


2059 


CD 

LD 

> 


O 

LU 

— D 

CD 

O 


03 

co 

> 

"O 

3 

"co 

C/3 


c/3 

03 

> 

1 

lo 

o 

Cl) 


o 

CO 

XJ 

cz 

X 

CO 

C/3 

CD 

o 

E 

00 

■ zz 

T3 

03 

C 

O 

3 

CO 

O 

CO 

CO 

_c 

£ 

o 

3 

o 

Q 

£ 

CO 

LU 

c 

cd 

"cd 

o 

>N 

“O 

Q- 

o 

"O 

=3 

“O 

03 

CO 

lo 

CD 

*c 

CD 

~CD 

03 

CL 

*o 

E 

C/3 

o 

O 

E 

03 

CO 

03 

Q_ 

X 

CO 

a) 

o 

£ 

C 

03 

SZ 

o 

o 

CL 

CO 

tr 

c 

Q 

o 

o 

LU 

a_ 

C/3 

O 

Q- 

3 

CO 

a> 

a. 

CO 

“O 

£ 

CD 

03 

c 

o 

co 

o 

CO 

X3 

03  „ 

$ 

CO 

c C 

03 

3 

=3  2 

c 

C/3 

11 

c 

3 

03 

Q 

> S 

■a 

LU 

c 

o 

X3  03 
03  "O 

$ 

03 

_c 

c3  c 

o 

03“ 

C/3 

03 

CD  CO 

Q 

CO 

C 3 

Cig> 

CO  c 

3= 

3 

CO 

03 

03 

> 

03 

CL  O 

£ 

1§ 

O ° 
*03  ^ 

“O 

*co 

X? 

^ ^ i 

03  1 

of 

3 

00 

"co 

> 

O 

o 

-£=s  x: 

- 

T— 

CM 

CO 

2060 


OBJECTIVES 


2061 
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ACCOMPLISHMENTS  IN  1 997 
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PROCEDURE  FOR 
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HSCT  ARROW  WING  M2.4-7A 
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2.  4%-Scale  M2.4-7 A/Tunnel  CFD  3.  4%-Scale  M2.4-7 A/Tunnel/Posts  CFD 
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OVERSET  GRID 

FOR  THE  4%-SCALE  HSCT  ARROW  WING  M2.4-7A 
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MAGGIE  CALCULATED  OVERSET  GRID 
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MAGGIE  CALCULATED  OVERSET  GRID 

FOR  THE  4%-SCALE  HSCT  ARROW  WING  M2.4-7A 


pMMBiBBBS 

!!!!■!!■■«  ygmB 

I IIIIIIIIIS 

i iiimniiHi  SSSniuiiM 

iiiiiiiBiBSiBiSpfSfi’;  ! ■ 

wmamtA 

iniiiiiiisMMratVhm 


1 3SS®8S5BftSS*S!!SSKgSSsSS8Sga™ilMIIIII 
I ^SSstesteiKJssiiBJMsaiBBiJiiiimniiii  i I 

l^iltWriMUSUIKHiMliJIBB 

l|i!jUi!U»!llii!ISS!IH:im8l!ill 
MsiMiWttKlHUiMfMIIII 


iiniiiiiai^i^r^iH 
llllllllll£»»P£Plf.!¥nl  .... 

i iiiiiiaaisiafifjfij'.ijf !«  i'smiii!  liiiiMSiui 

I lllllimiilWJlWM'IMll II  Kill  II.IKlIIII 
itiiiiiniafi»iRKrariKssir.Sji!iiiM!Siii!ii*s|i 
I IIIIIIIIIIPIfi2l'ill(S!'if  li!!S  iillll  IC.Sil  'llli 

iiiiiiiiiii£ariiis:!r.iiuimi!!.iiiH(!!irii; 
llllllllllimilllllliiiillll  liffiimM 

iiiiiiiiiiiiiiiiiiiimjimiMMIH 
iiiiiiiiiiimiiiiiiiiiiiiiii 
11111111111111111111111111!!! 


iiiiiiiiiin 

iiiiiiiiini 

llllllllllll 


iiiriiiiwimisuiiriiiKisiMiaaaaaaiiiii 1 

.illLt’il'IIIHItIMIIMISi'IBIJaflBBailllll  1 1 

riKiiiiiniiHMisiiiKiiMisaiaiiiiiiiiii' 

■illllMSIHIMAlllliKlIJJKifllMIIIIIIIII1 

i i;i  is  i uiiBihf  aa  iiiiui  mi  ii?aaiiaaaaaaaaiiin 
:ii<iii>!iiiuiKimi2»siiniRsiiiiii  mu  i 
iiikiiitiiiiik'iiiiiiiiiiiiiiiiiiimn  linn 
iiiiiiiiiiiiiiiiiiiimimiiiiiiami  inn  i 
iiiiiiiiiiiiiiiiiiiimiiiimiiii  hi  mii 
iiiiiiiiiiiiiiiiiimimmiiiiii  in  mu 
lllllllllllllllllll  II llllllllllll  in  mi 


2073 


CO 

LLI 

O 

< 

cc 

LU 


O 


< 

Q_ 

Q 

LL 

O 


CD 


T3 

CD 

a3 


W 

00 


CD 

E 

cc 


CD 

CT)-Q 

■E  5 
§:  cd 

X3 

"O  ci 

o 

F 7=5 

c <D 
CD  CO 
_d  CO 
CD 


CO 


CL 

CL 


_>•*  CO 

CO  q5 
JD  o 
O -O 
<0  CO 

o 

11 

£ x 

■p  CD 

II 

° CD 
CD  C 


JC  J_ 

I 2 

CD  ° 

o — ^ 
o 

4^  CO 
CD  •*- 


CD 


o O 

CD  "— 
JC  X 
CD 

o-g 

c:  > 


O) 


CO 

CD 


CT> 


CD 


CD 

g 

O 

CD 


"O 

c 

CO 

CD 

CL 

CL 


CD 


CO 

CD 

O 

TD 

CD 

C 

_o 

CO 

15 

5 

CD 


-o 

c 


2074 


CFD  PARTICLE  TRACES 

M2.4-7A  ARROW  WING 
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SURFACE  PRESSURE  CONTOURS 
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SURFACE  PRESSURE  CONTOURS 

ON  THE  UPPER  SURFACE  OF  ARROW  WING  M2.4-7A 
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Model/Tunnel/Posts 


SURFACE  PRESSURE  CONTOURS  ON  POSTS  AND  CROSS  SECTIONS 
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SURFACE  PRESSURE  CONTOURS  ON  MODEL/TUNNEL  WALL/POSTS 
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SURFACE  PRESSURE  CONTOURS 

TUNNEL  WALUPOSTS/4%-SCALE  M2.4-7A  MODEL 
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LOCATION  OF  THE  MOMENT  AXIS 
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LOCATION  OF  THE  MOMENT  AXIS 

INSIDE  NASA/AMES  12-FT  TUNNEL 
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COMPARISONS  OF  FORCES  AND  MOMENT 
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Pressure  range  1-  6 atmospheres. 

Supported  by  - MDC,  BCAG,  LMAS,  ARC,  LaRC. 
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mount  is  similar  in  diameter  to  the  post  mount  used  for  the  Ref  H 
Upflow  and  Interference  (U  & I)  test  in  the  14’x22’  tunnel,  and  the  U 
& I test  did  indicate  that  the  post  support  will  have  a significant  effect 
on  the  pitching  moment. 
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discrepancy  is  unclear. 

It  can  be  seen  from  the  experimental  pitching  moment  curves  that  the  shift  in 
pitching  moment  with  increasing  Reynolds  number  seen  in  the  tail-on  data 
presented  earlier  is  not  seen  here  with  the  tail  off. 
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Reynolds  number  variation.  It  is  therefore  clear  that  the  effects  of  the 
aeroelastic  model  deformation  should  be  fully  accounted  for  in  evaluating  the 
Reynolds  number  effects. 
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Surface  pressure  distribution  is  very  well  predicted  by 
CFL3D  Navier-Stokes  method 
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High  Lift  Technology  Development  (Task  33) 


Increase^  L/cy,  Qevejop  Analysi^Design  Methodology 


Goals 

Objectives 

Challenges 

Approaches 


Demonstrate  Greatly  Increased  L/D  Relative  to 
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The  powered  TCA-2  test  performed  at  Langley’s  14’  x 22’  resides  under  the  Test  Programs 
and  Techniques  Approach  within  Task  33,  High  Lift  Technology  Development. 
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Outline 


• Introduction  - LaRC  455  in  the  TCA  WTT  Plan 

• Test  Objectives 

• Test  Configuration  Matrix 

• 5%  TCA  Model  Description 

• Results 

• Conclusions 

• Lessons  Learned 

• Recommendations 
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The  test  was  performed  from  August  29th  to  the  9th  of  September,  1997.  The  average  run 
rate,  including  model  installation  was  1.5  runs  per  hour.  The  picture  above  depicts  the  aft 
end  of  the  powered  5%  HSCT  TCA  model.  The  model  will  be  discussed  in  more  detail  later 
in  the  presentation.  The  above  picture  will  be  referred  to  at  that  time. 
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WIND-TUNNEL  TEST  SCHEDU 

T C A MODEL 


; Ugh  Lift  Aerodynamics,  Long  Beach 

L E FOR  5 % 


19  9 7 

19  9 8 

JFM  AM  JJASON  D 

JFM  AM  JJASO  N D 

□ 

T C A - 1 
1 4 ’x  2 2 


T C A -2 
1 4 'x  2 2 


□ 


T C A -3 
Ames  12’ 


□ 

T C A -4 
14x22’ 


CZ3 


T C A - 5 
Ames  12’ 


The  TCA-2  wind-tunnel  test  was  the  second  in  a series  of  planned  tests  utilizing  the  5% 
Technology  Concept  Airplane  (TCA)  model.  Each  of  the  tests  was  planned  to  utilize  the 
unique  capabilities  of  the  NASA  Langley  14’x22’  and  the  NASA  Ames  12’  test  facilities,  in 
order  to  assess  specific  aspects  of  the  high  lift  and  stability  and  control  characteristics  of  the 
TCA  configuration.  However,  shortly  after  the  completion  of  the  TCA-1  test,  an  early 
projection  of  the  Technology  Configuration  (TC)  identified  the  need  for  several  significant 
changes  to  the  baseline  TCA  configuration.  These  changes  were  necessary  in  order  to  meet 
more  stringent  noise  certification  levels,  as  well  as,  to  provide  a means  to  control  dynamic 
structural  modes.  The  projected  changes  included  a change  to  the  outboard  wing  (increased 
aspect  ratio  and  lower  sweep)  and  a reconfiguration  of  the  longitudinal  control  surfaces  to 
include  a medium  size  canard  and  a reduced  horizontal  tail.  The  impact  of  these  proposed 
changes  did  not  affect  the  TCA-2  test,  because  it  was  specifically  planned  to  address  power 
effects  on  the  empennage  and  a smaller  horizontal  tail  was  in  the  plan  to  be  tested. 
However,  the  focus  of  future  tests  was  reevaluated  and  the  emphasis  was  shifted  away  from 
assessment  of  TCA  specific  configurations  to  a more  general  assessment  of  configurations 
that  encompass  the  projected  design  space  for  the  TC. 
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Test  Objectives  of  LaRC  455 


1.  Determine  Exhaust  effects  on  the  Empennage  in  Free-air 
(F-A)  and  in  Ground-Effect  (G.E.) 

• Baseline  (800  ft.2)  & Half  (400  ft.2)  Stabilizers 

• Takeoff  vs  Flow-thru  Eng.  Exhaust  Nozzles 

• Interference  Effects 

a.  High  Pressure  Air  Supply  Lines 

b.  Domed  Nacelles 


The  test’s  main  objective  was  to  quantify  the  effects  of  simulated  engine  plumes  on  the 
empennage  at  several  thrust  levels  with  emphasis  on  the  thrust  level  that  most  closely 
represents  take-off  conditions.  Additional  objectives  included  investigating:  1)  the  plumes’ 
effects  on  the  half-sized  (400  ft.2)  baseline  tail,  2)  difference  in  nozzle  positions  to  acquire 
the  desired  exhaust  effects,  and  3)  the  interference  effects  of  the  high  pressure  air  supply 
lines  and  domed  nacelles. 
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Test  Configuration  Matrix 


Five  configurations  that  were  tested: 

1.  Takeoff  flaps  (LE/TE  = 30/20),  3 Power  settings 

• Baseline  stabilizer 

• Half-size  stabilizer 

2.  Undeflected  flaps  (0/0),  3 Power  settings,  baseline  stabilizer 

3.  Takeoff  flaps,  flow-thru  nacelles 

• with  high  pressure  air  lines  installed 

• with  high  pressure  air  lines  removed 


The  test  configuration  matrix  was  structured  to  ensure  the  test’s  main  objective  would  be 
met  and  the  remaining  objectives  were  prioritized.  This  matrix  was  designed  into  the  plan  of 
test  (POT)  which  can  be  obtained  from  the  ‘ADAPT”  web  page.  The  High  Lift  and  Stability 
& Controls  objectives  for  best  trimmed  lift,  L/D,  simulation  and  control  power  assessment  of 
the  exhaust’s  effects  on  the  tail  were  combined  and  prioritized. 

The  test  was  centered  around  the  terminal  area  takeoff  and  climbout  conditions.  Most  of  the 
test  was  conducted  with  the  wing  in  the  takeoff  flap  configuration  - leading  edge  flaps  at  30  degrees 
and  trailing  edge  flaps  at  20  degrees  (LE/TE  30/20).  The  clean  wing  configuration  (0/0)  was  also 
studied.  The  undeflected  flap  configuration  was  tested  for  baseline  simulation  data  and 
additional  confirmation  about  downwash  at  the  empennage.  The  data  obtained  from  this  test 
will  be  additive  to  the  unpowered  TCA-1  test  data  to  culminate  in  a powered,  ground  effect 
and  free-air  TCA  database. 

The  test  included  runs  to  determine  interference  effects  that  must  be  accounted  for  because 
of  the  air  supply  lines  and  domed  inlets.  At  the  start  of  the  test,  the  need  to  test  the  ‘flow- 
thru’  nozzles  was  eliminated  by  the  ability  to  obtain  flow-thru  NPRs  with  the  takeoff 
nozzles. 

The  baseline  TCA  horizontal  tail  (HI)  was  primarily  used  to  determine  the  thrust  effects  on 
the  empennage.  HI  was  tested  with  both  0 degree  and  -30  degree  elevators.  A horizontal 
tail  with  half  the  area  of  HI  (H2)  was  also  tested.  H2  did  not  have  provisions  for  elevator 
deflections. 
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Tunnel  Test  Conditions 


Rec  = 7.54  x 106 
Mach  Number  of  0.245 
Dynamic  Pressure  of  85 1 
Nozzle  Pressure  Ratios 

h/b  (height*)  of  0.18  (14”),  0.21,  0.30,  0.49,  & 1.0  (77”) 
* height  (HGT)  is  W.L.  of  l.e.  50%  cbar  to  tunnel  floor. 


The  TCA-2  test  was  performed  at  the  LaRC  14’  x 22’  atmospheric  wind  tunnel.  The  TCA-1 
was  previously  tested  in  the  same  tunnel  with  a balance  located  in  both  the  body  and  the  tail 
section. 

Three  NPRs  were  tested  and  the  model’s  height  was  varied  as  indicated  above. 
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• Cover  Plate  and  Floor  Tufts  • Non-Metric/Metric  Bal.  Split 

• Hose  Slot  • HP  Hoses 

• Domed  Inlets  • Takeoff  Nozzles 


The  model’s  configuration  is  based  on  the  TCA  1080-1450  outer  mold  line  (OML) 
definition  released  in  March  of  1996.  The  model  has  been  designed  to  be  tested  in  the 
NASA  Langley’s  14’  by  22’  tunnel  both  powered  and  unpowered  and  also  in  NASA  Ames’ 
12’  pressure  tunnel,  unpowered,  for  Reynolds  number  effects.  The  model  was  constructed  to 
determine  the  TCA  configuration’s  low  speed  aerodynamic  performance,  provide  a database 
for  subsequent  analysis  and  simulation,  and  examine  the  simulated  nozzles’  propulsion 
effects  on  the  TCA  configuration’s  performance,  stability  and  control  in  ground  effect. 

The  model’s  wing  was  originally  designed  to  accommodate  high  pressure  air  plumbing 
necessary  for  powered  testing.  All  four  nacelles  were  supplied  with  separately  controlled 
high  pressure  air,  regulated  by  two  choke  plates  in  each  nacelle,  in  order  to  acquire  desired 
Nozzle  Pressure  Ratios  (NPR).  The  above  figures  depict  the  provisions  made  for  the  high 
pressure  air  supply  lines;  their  stiffness  introduced  concerns  about  tare  and  routing.  The  six- 
component  internal  balance  was  located  between  the  mid-body  and  empennage.  This 
facilitated  the  measurement  of  the  aft-body’s  control  surfaces’  aerodynamic  forces  without 
requiring  a complex  tare  reduction  mechanism  or  mathematical  tare  removal.  High  pressure 
air  supply  hoses  required  special  routing  through  a sliding  slot  in  the  floor  aft  of  the  support 
post.  A diagram  of  the  nacelle’s  called  a ‘Jet  Flow  Simulator’  follows  in  a few  pages. 

Two  sets  of  nozzles  were  built  for  the  nacelles.  One  set  was  intended  for  use  at  the  takeoff 
NPR  settings,  while  the  other  was  for  the  “flow-through”  setting.  Early  studies  showed  that 
swapping  the  nacelles  during  the  test  was  unnecessary  and  the  “flow-through”  NPR  setting 
could  be  run  with  the  takeoff  nozzles.  Slide  #4  can  be  viewed  for  additional  detail. 
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.2  static  pressure  probes 


5 total  pressure  probes 


r 


2 thermocouples 


Nozzle  exit  area 
= 8.76  sq.  in. 


> Constant  Area  Rectangular  Section  = 9.83  sq.  in.  < 


Jet  Flow  Simulator 


The  Jet  Flow  Simulator  diagram  depicts  its  functional  construction.  The  high  pressure 
supply  lines  injected  air  at  approximately  340  psi  to  the  nacelle’s  plenum.  The  pressurized 
air  was  moved  along  the  circular  to  rectangular  transition  and  choked  through  the  two  choke 
plates  #1  and  #2  to  achieve  the  desired  internal  flow  conditions.  The  Nozzle  Pressure  Ratio 
was  determined  by  dividing  the  total  pressure  at  the  charging  station  over  the  ambient  static 
pressure.  This  enabled  an  exhaust  plume  that  most  closely  represented  the  full-scale  plume. 

Mass  flow  rate  was  determined  by  pressures  and  temperatures  obtained  from  the  total  and 
static  pressure  probes  and  the  thermocouples.  Total  pressure  probe  placement  was 
important.  Erroneous  total  pressure  readings  could  have  been  obtained  if  the  estimated 
internal  flow  profile  was  not  known  prior  to  start  of  the  test.  The  probes  were  placed  in  a 
cosine  distribution  to  best  capture  the  conventional  channel  flow  profile. 

The  takeoff  and  flow-thru  nozzles  were  separate  pieces  attached  to  the  aft  end  of  the 
nacelles.  The  designed  nozzle  deflections  were  similar.  Since  the  flow-thru  NPR  was 
attainable,  the  time  spent  to  change  between  the  two  nozzles  was  saved  by  deleting  the  flow- 
thru  nozzle  from  the  test  matrix. 


af_revu  2/25/98 


2170 


Airframe  Review  2/98 


High  Lift  Aerodynamics,  Long  Beach 

Results 


ACm  Exhaust_gefa  (Cm  full  thrust 

* (Cm  full  thrust 


m flow 


-thru) 

hru) 


ground  effect 


m flow-thru/  free  air 


“If  ACn,  Exhaust_gefa  > 0.015,  test  entire  matrix.” 


1 . Ground  Effects  achieved  by  varying  the  model  HGT  from 
77”-*-  14”. 

2.  Exhaust  Effects  via  varying  NPR:  1 .04  & 1.51. 

The  following  data  are  presented  as  plots  of  pitching  moment  coefficient  both  absolute  and 
incremental.  Pitching  moment  was  selected  as  the  parameter  to  be  used  as  a metric  because  small 
changes  in  the  tail’s  local  flow  field  are  easily  observed.  Although  the  results  are  shown  and 
discussed  as  incremental  changes  in  pitching  moment,  the  results  are  also  quantified  in  terms  of  the 
change  in  the  local  angle  of  attack  at  the  tail  resulting  from  a change  in  either  stabilizer  incidence, 
downwash,  or  airplane  angle  of  attack.  The  change  in  local  angle  of  attack  at  the  tail  was 
determined  by  dividing  the  measured  incremental  pitching  moment  by  the  stabilizer  effectiveness 
(typically  -0.0045/deg.).  Although  the  above  equation  sums  up  the  additive  effects  of  both  ground 
effect  and  power,  the  following  plots  predominantly  separate  these  effects  for  clarity.  The  plots  and 
related  text  predominantly  address  the  10  degree  angle  of  attack  region.  This  was  chosen  as  a 
reasonable  alpha  for  takeoff  and  initial  climbout  regimes. 

A delta  of  0.015  in  pitching  moment  coefficient  was  agreed  on  by  team  researchers  prior  to  the  start 
of  the  test  in  order  to  determine  how  much  of  the  test  matrix  would  have  to  be  tested.  If  the  delta 
obtained  was  greater  than  0.015,  the  test  was  to  progress  through  the  entire  test  matrix  of  heights 
and  power  settings.  The  delta  obtained  early  on  in  the  test  was  greater  than  the  0.015;  therefore  the 
test  proceeded  accordingly. 

Testing  in  ground  effect  was  performed  by  moving  the  model  from  the  tunnel  centerline  of  77”  to 
14”  above  the  tunnel  floor.  An  NPR  of  1 .26  was  also  obtained  but  used  mainly  for  trends  and  not 
presented  here. 

The  test  centered  around  the  takeoff  flap  setting  of  LE/TE  = 30/20  with  the  stabilizer/elevator 
deflected  to  -15/-30.  This  stabilizer/elevator  deflection  provides  the  most  airplane  nose-up  pitching 
moment  simulating  the  nosewheel  lift-off  maneuver.  Although  the  data  was  obtained  on  the  aft 
balance,  all  pitching  and  yawing  moments  are  referenced  to  the  50%  cbar  of  the  wing.  Tail-off  and 
nominal  stabilizer  settings  of  -15/0  were  also  tested  to  determine  stabilizer  effectiveness  and 
downwash  at  the  horizontal  tail. 
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Pitching  Moment  Variation  with  Angle  of  Attack 
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The  above  plot  shows  the  impact  of  domed  inlets  at  both  77  inches  and  14  inches  ground 
heights.  A positive  increment  in  pitching  moment  due  to  the  presence  of  the  domes  is  seen 
at  both  heights.  While  the  impact  of  the  domes  is  relatively  small  at  the  maximum  height, 
equivalent  to  a 1/4  degree  change  in  tail  angle  of  attack,  the  impact  at  the  minimum  height  is 
significantly  more,  about  3/4  degree  change  in  tail  angle  of  attack. 
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Pitching  Moment  Variation  with  Angle  of  Attack 
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The  above  chart  shows  the  impact  of  the  air  supply  lines  on  the  measured  pitching  moment 
coefficients.  A positive  increment  in  pitching  moment  is  seen  due  to  the  presence  of  the  air 
supply  lines  in  free-air,  while  a negative  increment  is  seen  in  ground  effect.  The  magnitude 
of  the  increments  at  10  degrees  angle  of  attack  appear  to  be  the  same  at  both  heights  and  are 
equivalent  to  a 1/4  degree  change  in  tail  angle  of  attack. 
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Angle  of  Attack 


The  above  chart  shows  the  impact  of  NPR  on  pitching  moment  in  both  free-air  and  ground 
effect  for  the  takeoff  flaps  and  the  stabilizer/elevator  deflection  of -15/-30  configuration.  As 
NPR  is  increased  from  1.04  to  1.51,  a significant  increase  in  pitching  moment  is  observed  at 
the  minimum  height,  approximately  equivalent  to  a 3.5  degree  change  in  stabilizer 
incidence.  However,  at  the  maximum  height,  the  effect  of  NPR  is  greatly  diminished  and 
results  in  a decrease  in  pitching  moment,  approximately  equivalent  to  a 3/4  degree  change  in 
stabilizer  incidence. 


The  overall  nose-down  tendency  in  free-air  is  most  likely  due  to  the  entrainment  of  the 
wing’s  downwash  with  the  exhaust  plume  thus  reducing  the  local  angle  of  attack  of  the 
horizontal  tail.  The  crossover  in  free-air  at  16  to  18  degrees  angle  of  attack  indicates  that 
increasing  NPR  has  little  impact  on  pitching  moment  because  the  tail  is  in  the  wing’s  wake. 

The  ground  effect  set  of  curves  show  a significant  amount  of  nose-up  due  to  a dynamic 
pressure  or  venturi  effect  acting  on  the  lower  surfaces  of  the  horizontal  tail.  The  slide  on 
page  19  shows  oil  flow  visualization  highlighting  this  phenomenon. 

The  next  slide  summarizes  the  impact  of  the  interference  effects  on  the  above  measured  data. 
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The  above  chart  shows  the  combined  effects  of  domed  inlets  and  air  supply  lines  on  the 
measured  incremental  pitching  moment  coefficients.  The  final  increment  data  (not  including 
wind  tunnel  data  reduction  schemes)  will  be  obtained  by  removing  the  cumulative  power 
artifact  effects  (interference  of  domed  inlets  and  air  supply  lines)  from  the  measured  data. 
The  final  power  effect  increments  in  both  free-air  and  ground  effect  can  be  applied  to  the 
simulation  database. 

Note  that  in  ground  effect  the  artifact  effects  negate  each  other  while  in  free-air  they  are 
additive  and  thus  increase  the  measured  pitching  moment  increments.  Final  incremental 
data  for  the  above  case  shows  the  free-air  pitching  moment  increment  decreasing  to  just  over 
one  degree  of  tail  angle  of  attack  while  the  ground  effect  increment  remains  the  same. 
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Pictures:  Top  1 ->  2,  Bottom  3 ->  4. 

The  above  pictures  illustrate  some  of  the  oil  flow  visualization  performed  during  the  test. 
These  pictures  are  of  the  same  run  from  different  perspectives.  Note  that  the  oil  pattern  in 
picture  2 looks  very  symmetrical  but  when  viewed  from  a different  angle,  such  as  picture  1 , 
the  flow  lines  are  somewhat  distorted  and  seem  ‘braided.’  Although  it  was  thought  that  the 
air  lines  were  going  to  cause  significant  interference,  they  did  not  by  virtue  of  the  results 
shown  on  the  previous  chart.  Pictures  1-3  depict  the  wing’s  cross-flow  and  vortices  even  at 
this  relatively  small  angle  of  attack  of  8 degrees.  It  should  be  evident  from  the  above 
pictures  that  the  air  is  sped  up  under  the  tail  consistent  with  a venturi  effect  which  creates 
lower  pressures  and  thus  results  in  the  increased  airplane  nose-up  pitching  moment 
increment.  Picture  4 is  just  a side  view;  the  tail’s  upper  surface  can  be  seen. 
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Pitching  Moment  Variation  with  Angle  of  Attack 


Angle  of  Attack,  a 
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The  next  four  slides  show  the  effect  of  NPR  on  both  elevator  and  stabilizer  effectiveness  for 
stabilizer/elevator  deflections  of  -15/0  and  -1 5/-30  with  takeoff  flaps  LE/TE  30/20.  The  data 
are  presented  in  free-air  and  in  ground  effect  separately.  The  data  are  presented  at  the  two 
heights  as  both  absolute  and  incremental. 

The  above  chart  shows  the  absolute  pitching  moment  data  in  free-air.  As  NPR  is  increased, 
the  elevator’s  effectiveness  is  shown  to  reduce  by  12%  while  the  stabilizer’s  effectiveness  is 
increased  by  4.5%.  This  reduction  of  elevator  effectiveness  is  puzzling  and  is  better 
illustrated  in  the  next  slide. 
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Incremental  Pitching  Moment  Variation  with  Angle  of  Attack 


The  above  chart  shows  the  incremental  pitching  moment  effect  of  increasing  NPR  in  free-air 
on  the  two  specific  stabilizer/elevator  combinations.  The  increment  is  calculated  by 
subtracting  the  absolute  data  at  NPR  = 1.04  from  the  absolute  data  at  1.51  for  each 
stabilizer/elevator  configuration  shown.  The  incremental  pitching  moment  effect  of  NPR  on 
the  -15/0  stabilizer  deflection  is  positive  over  the  entire  angle  of  attack  range.  The  positive 
increment  equates  to  approximately  a -2  degree  stabilizer  incidence  change. 

Note,  as  NPR  is  increased  with  the  stabilizer/elevator  deflection  of -15/-30.  the  pitching 
moment  increment  is  negative  for  the  terminal  area  alpha  range.  This  reduction  in  elevator 
effectiveness  equates  to  approximately  a 3/4  degree  change  in  stabilizer  incidence. 

Repeat  runs  are  also  shown  in  this  plot.  Repeatability  is  shown  to  be  reasonable,  particularly 
at  1 0 degrees  alpha.  Consistent  NPR  settings  would  obviate  much  of  the  variance  in  the  data 
and  will  be  discussed  later. 


afjnevu  2/25/98 


2178 


Airframe  Review  2/98 


Hfrjh  Lift  Aerodynamics,  Long  Beach 


Pitching  Moment  Variation  with  Angle  of  Attack 


Angle  of  Attack,  a 
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The  above  chart  shows  the  absolute  pitching  moment  data  in  ground  effect.  This  chart 
shows  that  increasing  NPR  in  ground  effect  has  a beneficial  effect  on  both  stabilizer  and 
elevator  effectiveness.  This  significant  increase  in  stabilizer  effectiveness  equates  to 
approximately  a -4  degrees  change  in  additional  stabilizer  incidence.  The  increase  in 
elevator  effectiveness  is  similar  to  a -2  degree  change  in  stabilizer  incidence.  The  following 
slide  better  illustrates  the  benefits. 
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The  above  chart  shows  that  increased  NPR  increases  both  the  stabilizer  and  elevator 
effectiveness  in  ground  effect  for  all  positive  angles  of  attack.  The  increments  are  calculated 
by  subtracting  the  absolute  data  at  NPR  = 1.04  from  the  absolute  data  at  1.51  for  each 
stabilizer/elevator  configuration  shown. 

Repeat  runs  are  also  shown  in  this  plot.  Repeatability  of  the  data  is  worse  in  ground  effect 
than  in  free-air  but  is  sill  acceptable,  particularly  for  the  -15/-30  stabilizer/elevator 
configuration.  Overall,  repeatability  could  have  been  improved  by  maintaining  more 
consistent  NPRs.  The  tail  balance  showed  no  significant  anomalies  during  the  test  that 
would  have  resulted  in  unsatisfactory  repeatability. 
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Incremental  Pitching  Moment  Variation  with  Angle  of  Attack 


Angle  of  Attack,  a 
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The  above  chart  shows  the  incremental  effect  of  NPR  on  the  half-sized  (H2)  horizontal  tail 
(400  sq.  ft.)  both  in  free-air  and  ground  effect.  The  incremental  effect  of  NPR  on  the 
baseline  (HI,  800  sq.  ft.)  is  also  shown  for  comparison.  Both  HI  and  H2  are  at  a 
stabilizer/elevator  deflection  of  -15/0  in  the  above  incremental  data  plot. 

In  general,  the  incremental  effect  of  NPR  on  the  H2  tail  in  free-air  is  more  than  half  the 
increment  observed  on  the  HI  tail  at  low  to  moderate  angles  of  attack.  This  incremental 
effect  is  approximately  half  at  higher  angles  of  attack  until  20  degrees  is  reached  where  the 
effect  reduces  to  no  observable  effect  at  26  degrees  alpha. 

However,  in  ground  effect  the  incremental  effect  of  NPR  on  the  H2  tail  is  approximately  half 
that  observed  for  the  HI  tail  up  to  about  4 degrees  angle  of  attack.  Beyond  4 degrees  alpha, 
the  incremental  effect  on  the  H2  tail  begins  to  rapidly  diminish  with  no  effect  evident  at  1 0 
degrees  alpha. 
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Angle  of  Attack,  a 


The  above  plot  shows  the  effect  of  NPR  on  the  0/0  and  30/20  flap  configurations  at  both 
free-air  and  ground  effect  heights  for  the  incremental  pitching  moment  of  the 
stabilizer/elevator  deflected  -15/-30.  The  0/0  flap  deflection  was  tested  in  conjunction  with 
the  30/20  flap  deflection  in  order  to  provide  an  increment  for  the  simulation  database.  The 
reduced  level  of  pitching  moment  for  the  0/0  flap  when  compared  to  the  30/20  flap 
configuration  is  the  result  of  reduced  downwash  and  thus  reduced  local  alpha  at  the  tail. 

The  trends  in  pitching  moment  increment  due  to  NPR  and  height  changes  are  similar  for  the 
0/0  flaps  when  compared  to  the  30/20  flaps.  One  exception  is  when  NPR  is  increased  in 
free-air  for  the  30/20  flap  deflection;  -4  to  16  degrees  alpha  shows  an  airplane  nose-down 
pitching  increment  while  the  0/0  flaps  shows  the  opposite  increment.  The  ground  effect 
curves  of  the  two  different  flaps  settings  show  comparable,  positive  increments. 
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Y awing  Moment  Variation  with  Angle  of  Sideslip 
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Angle  of  Sideslip,  p 


The  above  plot  shows  the  effect  of  NPR  on  directional  stability  for  alpha  = 10  degrees.  Both 
in  free-air  and  in  ground  effect,  the  data  show  that  increasing  NPR  increases  directional 
stability.  While  the  effect  of  NPR  is  of  primary  interest,  ground  effect  is  shown  to 
significantly  reduce  directional  stability  particularly  at  small  sideslip  angles.  At  all  angles  of 
sideslip,  the  yawing  moment  due  to  sideslip  is  greater  in  free-air  than  in  ground  effect. 
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Conclusions 


• Impact  of  exhaust  on  tail  is  a ‘non-issue’  as  the  stabilizer  is  far 
from  stalled. 

• ACm_NPR  (30/20  Flaps):  negative  in  free-air  due  to  reduced 
downwash  from  plume  entrainment  and  positive  in  ground  effect 
due  to  a dynamic  pressure  increase  on  tail  near  the  ground. 

• Elevator  effectiveness  lessens  in  free-air  and  increases  in  ground 
effect  while  stabilizer  effectiveness  increases  at  both  heights. 

• Half-sized  horizontal  tail  effectiveness  due  to  NPR  was  greater 
than  half  of  the  baseline  tail’s  in  free-air  and  less  than  half  in 
ground  effect. 

• Increased  Cnp  due  to  NPR. 

The  main  conclusion  drawn  from  the  data  thus  far  indicates  that  the  exhaust  plume 
simulating  takeoff  power  shows  no  detrimental  effects  to  longitudinal  control  power. 

An  additional  nose-up  pitching  moment  increment  due  to  NPR  was  obtained  in  ground  effect 
for  the  -15/-30  stabilizer/elevator  deflection  for  both  flap  configurations  tested.  However, 
the  flaps  30/20  free-air  case  showed  a slight  airplane  nose-down  increment  with  the 
application  of  power. 

Stabilizer  effectiveness  was  shown  to  increase  in  both  free-air  and  ground  effect  for  the 
baseline  horizontal  tail  with  the  application  of  power.  Elevator  effectiveness  on  the  baseline 
tail  was  reduced  in  free-air  while  in  ground  effect  it  was  shown  to  increase.  The  half-sized 
horizontal  tail  effectiveness  increase  due  to  NPR  was  shown  to  be  greater  than  half  of  the 
baseline  tail’s.  This  effect  due  to  NPR  was  reduced  in  ground  effect.  Both  increases  in 
stabilizer  and  elevator  effectiveness  benefit  the  nosewheel  lift-off  maneuver  at  takeoff. 

Power  increased  directional  stability. 
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Lessons  Learned 


• High  pressure  flex-hoses  don’t  have  to  be  stowed. 

• Repeatability  from  precise  control  of  NPR  a must. 

a.  NPR  control  of  0.05  desired  while  0.1  was  the  norm. 

b.  NPR  settings  would  drift  as  run  was  on-going. 

c.  No  acceptable  repeatability  by  human  intervention. 

• Leak  rate  less  than  0.005  lbs/sec  difficult  to  meet. 

• Precise,  consistent  stabilizer  drive  needed. 

• Good  mass  flow  calculation  is  f(Probe  placement  within  the 
nacelle). 

• Prominently  mark  pertinent  model  pieces. 

• Probes  fail. 

Prior  to  the  start  of  the  test,  the  air  supply  lines  were  anticipated  to  have  a significant  effect 
thus  reducing  the  quality  of  the  data.  This  turned  out  not  to  be  the  case  as  shown  previously. 

Data  repeatability  could  be  improved  by  maintaining  more  consistent  NPR  settings.  NPR 
settings  varied  during  runs  and  were  not  easily  maintained  even  with  the  assistance  of  human 
intervention.  Better  control  of  NPR  would  reduce  any  potential  repeatability  problems. 

Maintaining  a leak  rate  below  the  target  of  0.005  lbs/sec  was  difficult.  The  model  was  very 
susceptible  to  leaking  if  touched  by  the  mechanics. 

The  stabilizer  was  susceptible  to  drift.  Maintaining  constant  stabilizer  incidence  would 
improve  repeatability. 

The  placement  of  total  pressure  probes  affect  the  mass  flow  calculations.  The  flow  profile 
must  be  known  apriori  in  order  to  place  total  pressure  and  temperature  probes  in  a logical  set 
of  locations  for  good  mass  flow  calculations. 

A portion  of  the  test  was  performed  with  the  wrong  outboard  flaps,  which  cost  valuable 
research  time.  Model  parts  need  to  be  prominently  marked  for  easy  identification. 

Pressure  and  temperature  probes  failed.  This  led  researchers  to  substitute  readings  from  one 
probe  to  that  of  another  location.  This  leads  to  errors  in  specific  flow  quantities  and  can 
provide  misleading  results  when  data  are  reduced  for  mass  flow  calculations. 
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Recommendations 


• Perform  an  isolated  check  out  of  powered  nacelle(s). 

• Design  larger  powered  models  to  further  reduce  any  high  pressure 
air-line  routing  impact  if  the  testing  facilities  are  capable. 

• Ensure  tail  drive  mechanism  has  no  hysteresis  in  its  inherent 
design  by  increasing  estimated  loads. 

• Revamp  the  controllability  of  sustaining  NPR. 

• Acquire  ability  to  seed  plume. 

• Heated  plumes  recreate  buoyancy  (ejectors/bypass). 

• Know  that  ‘leaks’  are  not  a ‘good’  thing. 

The  following  recommendations  are  given  for  future  powered  tests.  It  seems  logical  to 
ensure  that  the  workings  of  the  jet  flow  simulators  (JFS)  are  fully  functional  before  the 
actual  test  start  date.  Perhaps  an  isolated  test  of  the  JFSs  is  necessary.  Leak  rates  could  have 
been  reduced  earlier  which  would  have  allowed  the  actual  data  taking  to  have  started  earlier. 

Data  repeatability  can  always  be  improved  by  reducing  the  mast  wake  turbulence,  increasing 
the  controllability  of  NPR  and  eliminating  the  drift  in  tail  incidence.  A fairing  around  the  air 
supply  lines  could  have  reduced  further  any  interference  from  these  artifacts.  Controllers 
used  to  monitor  NPR  need  to  be  analyzed.  As  mentioned  before,  NPR  settings  varied  greatly 
during  a test  run. 

Hysteresis  in  any  control  surface  drive  mechanism  needs  to  be  determined  by  calibrating 
with  substantially  more  loads  than  anticipated  in  the  test. 

There  is  a definite  need  to  view  the  exhaust  plume  in  three-dimensions.  Oil  flow  does  not 
give  height  nor  breadth.  This  would  have  been  more  insightful  and  should  be  provided 
where  flow  visualization  is  needed  or  requested. 

Heated  plumes  create  buoyancy  which  better  simulate  the  real  exhaust  plume.  The  effect  of 
ejector  and  bypass  door  turbulence  was  not  investigated  in  this  test.  These  are  added 
complications  but  factors  that  lead  to  a closer  approximation  of  the  exhaust  plumes.  Some 
of  these  items  were  looked  at  in  the  HEAT  1 test. 

In  closing,  when  performing  a powered  test,  know  that  leaks  in  pressurized  nacelles  do  not 
assist  in  gathering  good  data. 
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Aerodynamics  Applications,”  Bryan  Campbell. 

Dynamic  acoustic  detection  is  a technique  that  has  been  discussed  for  some  time 
now  but  has  not  yet  been  extensively  developed.  “Dynamic  Acoustic  Detection  of 
Boundary  Layer  Transition,”  by  Jonathan  (J.r)  Grohs  and  Guy  Kemmerly  discusses 
this  technique. 
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Dynamic  Acoustic  Detection 
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using  tms  tecnmque  would  be  limited. 

TSP  Imaging 

- A significant  amount  of  development  work  had  already 
been  completed  for  this  technique  in  academia. 

- The  technique  had  not  been  tested  in  a production  facil 
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addition,  the  length  of  pressure  tubing  had  been  found  to 
affect  the  data  acquired.  A substantial  effort  would  be 
required  to  determine  if  this  effect  could  be  overcome 
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Temperature  Sensitive  Paint  (TSP) 
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sensitivity  of  ~1%  change  in  luminescence  per  °F 
yielded  good  results.  Paints  with  a sensitivity  of 
-0.85%  change  in  luminescence  per  °F  did  not  yield 
consistent  results. 
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- Rincon  Mountain  Design  (not  currently  use 
• Redundant  image  processing  systems  exist: 

- PC  based  system  for  ratioing  images 

- SGI  based  system  for  complete  image  processing 
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Comparison  of  Sublimating  Chemical 
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No  evidence  of  lamina 
at  high-lift  design  poin 
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Mach  = 0.9 
Alpha  = +3  Degrees 
Upper  Surface 
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and  spotting  the  corresponding  laminar  fraction  point  at  that  drag  level.  This  would 
shift  the  analytical  skin-friction  curves  up  or  down  to  allow  extrapolation  without 
having  to  acquire  high  Reynolds  data. 

This  chart  also  illustrates  the  larger  than  expected  laminar  surface  area  present. 
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All  experimental  data 
represents  free  transition 
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All  experimental  data 
represents  free  transition 
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baseline,  conventional  trip  indicating  that  the  transition 
strip  was  undersized  for  some  of  the  conditions  tested. 
A trip  with  increased  grit  size  was  tested  which 
eliminated  this  clouding. 
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flow  was  observed  on  the  leading  edge  of  the  transonic 
configuration.  The  boundary-layer  is  tripped  at  the 
leading-edge  flap  hingeline. 
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This  presentation  includes  a summary  of  a recent  experimental  study  of  the 
static  and  dynamic  ground  effects  for  low  aspect  ratio  wings.  The  authors 
would  like  to  thank  the  many  members  of  the  Dynamic  Ground  Effects  (DGE) 
Team  whose  contributions  were  invaluable  in  this  effort. 
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High  Lift  Technology  Development  (Task  33) 


Goals 


Objectives 


Challenges 


Approaches 


Demonstrate  Greatly  Increased  L/D  Relative  to 
SST  Technology  (Suction  Parameter  £92%) 


Define  Preferred 
High  Lift  System 
~H ----- 


Technology 
Readiness  Level  £ 6 


Analysis/Design 

Methodology 


Aggressive 
Tech  Projection 


Validation 
(Little  Data) 


Concept 


Test  Programs 


Planfomri,  Viscous,  Scale, 
Ground,  Power,  Canard  Effects 


Program 


..... — t. 

_ Boundary  Layer 
Control 

■ 



-(Vortex  Raps  | 

- Powered  Testing  | 

]/Musui  iuu  r tuw  ridpb  | 

■ Dynamic  GE  Testing! 

H Programmed  Raps  | 

_ Pressure  & Temp.  | 

Analytical 

Methods 


HTechnology  Projection  | 


Sensitive  Paints 


jLinear  Methods  I 


•{  Navier-Stokes  Methods 


This  study  supports  the  Dynamic  GE  Testing  Program  under  Test  Programs 
and  Techniques  of  the  HSR  High  Lift  Technology  Development  Task. 
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HSCT  High  Lift  Aerodynamics 

Outline 

• General  Test  Information  (Owens) 

• Static  Ground  Effects  (SGE)  Data  Analysis  and 
Results  (Powell) 

• Dynamic  Ground  Effects  (DGE)  Data  Analysis 
and  Results  (Owens) 

• Comparison  of  Wind  Tunnel  and  Flight  Ground 
Effects  Data  (Curry) 


This  presentation  is  divided  into  four  main  sections.  First,  Lewis  Owens 
(Langley  Research  Center)  will  cover  some  general  Langley  14-by-22-Foot 
Subsonic  Tunnel  (14x22  ft)  test  information  (LaRC  Test  462  — October  ‘97). 
Next,  Art  Powell  (Boeing-Long  Beach)  will  present  some  of  the  analysis  and 
results  for  the  static  ground  effects  (SGE)  data,  which  provided  a baseline  for 
comparison  with  the  dynamic  ground  effects  (DGE)  data.  Lewis  will  follow 
with  a presentation  of  the  analysis  and  results  for  the  dynamic  ground  effects 
data.  Finally,  Bob  Curry  (Dry den  Flight  Research  Center)  will  present  a 
summary  and  status  of  recent  Tu-144  flight  test  ground  effects  results.  He  will 
also  compare  these  results  with  data  obtained  from  the  wind  tunnel  test. 
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Previous  ground  effects  data  (University  of  Kansas)  illustrates  the  potential  for 
a difference  between  DGE  and  SGE.  This  plot  of  the  percent  increase  in  lift  as 
a function  of  the  aspect  ratio  is  used  to  show  that  some  wing  planforms  exhibit 
significant  differences  between  DGE  and  SGE.  At  lower  aspect  ratios,  the 
DGE  lift  increase  may  only  be  50%  of  the  SGE  lift  increase.  If  this  happens  to 
be  the  case  for  the  HSCT  planforms  (aspect  ratio  currently  about  2),  then  the 
expected  SGE  lift  increments  may  over  predict  the  flight  (dynamic)  ground 
effect.  The  significance  of  this  over  estimation  is  that  the  flight  control 
surfaces  may  be  under  designed. 


University  of  Kansas  data  reference: 

R.  C.  Chang  and  V.  U.  Muirhead,  “Effect  of  Sink  Rate  on  Ground  Effect  of 
Low-Aspect-Ratio  Wings,”  Journal  of  Aircraft,  Vol.  24,  No.  3,  March  1987, 
pp.  176-180. 
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Test  Objectives 

• Shakedown  the  new  14’x22’  DGE  cart  and 
instrumentation  systems,  and  develop  the  overall 
system  as  a way  of  acquiring  DGE  data 

• Determine  the  extent  of  DGE  on  HSCT  planforms 

• Understand  how  planform  variables  such  as  aspect 
ratio  and  sweep  affect  DGE 


DGE  cart  hardware  existed  but  had  not  been  scheduled  for  use  in  the  14’x22’ 
tunnel.  The  HSCT  DGE/SGE  difference  question  provided  an  opportunity  to 
use  this  new  cart.  The  HSR  test  was  scheduled  and  a test  plan  was  developed 
to  meet  three  main  test  objectives.  First,  we  had  to  develop  experience  to  be 
able  to  effectively  use  the  DGE  cart.  Each  step  toward  operating  the  cart  in  the 
tunnel  (for  the  first  time)  involved  a major  “debugging”  effort.  This 
shakedown  process  would  also  include  validation  of  the  DGE  test  technique  by 
making  comparisons  with  the  Tu-144  flight  test  ground  effects  database. 
Second,  we  wanted  to  obtain  DGE/SGE  data  on  as  many  HSCT  planforms  as 
possible  to  try  and  provide  an  answer  to  the  DGE/SGE  difference  question. 
Finally,  we  wanted  to  test  a number  of  different  low  aspect  ratio  planforms  to 
be  able  to  understand  the  dominant  geometry  factors  that  may  contribute  to  the 
potential  DGE/SGE  difference. 
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Here  is  a side  view  sketch  of  the  DGE  cart.  The  basic  cart  is  similar  to  other 
14’x22’  model  carts.  The  differences  are  primarily  in  the  model  support 
structure.  A large  support  strut  is  hydraulically  controlled  (vertical  drive)  to 
vary  the  model  height  above  the  cart  floor.  Also,  the  strut  has  a hydraulic 
pitch  drive  that  makes  it  possible  to  change  the  model  attitude  during  a 
dynamic  plunge.  Finally,  although  not  used  in  this  test,  a yaw  drive  allows  the 
model  to  be  yawed  with  respect  to  the  oncoming  flow  (prior  to  DGE  plunge) 
so  that  ground  effects  data  can  be  obtained  with  angle  of  sideslip. 
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Basic  Cart  Performance 

• Executes  pre-programmed  model  trajectory 

• Sink  rates  varying  from  0 to  15  ft/sec 

• Cart  floor  boundary  layer  minimized  with  tunnel 
boundary-layer  removal  system 

• Height  limits  are  89  in.  down  to  about  5 in. 
(depends  on  model  pitch) 

• Pitch  limits  are  -10  degrees  to  +50  degrees 
(depends  on  model  height) 


Here  are  some  of  the  basic  performance  parameters  for  the  new  DGE  cart 
design.  Note  that  we  did  have  some  operational  problems  with  the  cart. 
During  the  test,  we  were  only  able  to  reach  sustained  sink  rates  of  about  9 ft/ 
sec  due  to  temporary  “fixes”  made  to  the  hydraulic  control  system. 
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Models  Tested 

• Model  #10:  Elliptic  wing  planform  with  centerbody 
(AR=7.0)-  No  flaps  deflected 

• Model  #7:  Tu-144  wing  planform 
(AR=  1.627)-  TE  flaps  deflected  10  deg 

• Model  #6:  TCA  wing  planform 
(AR=2.027)  - TE  flaps  deflected  10  deg 


Originally,  the  test  plan  included  10  different  flat-plate  wing  planforms  (one 
high  AR  wing  and  nine  low  AR  wings).  Note  that  none  of  the  wing  planforms 
had  wing  twist  or  camber.  The  operational/training  problems  experienced 
with  the  new  cart  contributed  to  the  reduction  in  the  number  of  models  that 
were  tested.  Only  three  models  were  tested.  Model  #10  was  used  to  get  an 
indication  of  the  tunnel  flow  angularity  with  different  model  heights.  This  was 
done  both  statically  and  dynamically.  Model  #7  provided  an  opportunity  to 
validate  the  wind  tunnel  DGE  data  by  making  comparisons  to  actual  Tu-144 
flight  data.  Model  #6  provided  an  opportunity  to  assess  the  DGE/SGE 
differences  for  a current  HSCT  baseline  wing  planform. 
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HSCT  High  Ltft  A+rodyrmmtcs 

Elliptic  Wing  Planform 


This  is  a sketch  of  model  10  wing  planform  with  centerbody.  The  model  was 
tested  in  the  “high-wing”  configuration  such  that  the  centerbody  was  under  the 
wing.  The  dimensions  are  included  to  give  a better  sense  of  the  model  scale. 
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Tu- 144  Wing  Planform 


This  is  a sketch  of  the  model  7 Tu-144  wing  planform.  The  model  was  tested 
in  the  “low-wing”  configuration  with  the  balance/balance  fairing  mounting  to 
the  top  of  the  wing.  The  dimensions  are  included  to  give  a better  sense  of  the 
model  scale  and  location  of  trailing-edge  flaps. 
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This  is  a sketch  of  the  model  6 Technology  Concept  Airplane  (TCA)  wing 
planform.  The  model  was  also  tested  in  the  “low-wing”  configuration  with  the 
balance/balance  fairing  mounting  to  the  top  of  the  wing.  The  dimensions  are 
included  to  give  a better  sense  of  the  model  scale  and  location  of  trailing-edge 
flaps. 
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Wind  Tunnel  Measurements 


• Six  component  force/moment  balance 

• Six  model  accelerometers 

-3  linear  accelerations 
-3  angular  accelerations 

• Four  sting  accelerometers  (data  not  used) 

• Tunnel  Optotrak  system 

-model  height 
-model  pitch 


The  basic  measurements  provided  during  the  wind  tunnel  test  included  a force/ 
moment  balance,  accelerometers  and  an  optical  tracking  system.  The  six 
model  acceleration  measurements  were  used  in  combinations  to  provide  3 
linear  and  3 angular  accelerations  of  the  model  reference  center.  The  four 
sting  mounted  accelerometer  measurements  tracked  the  model  accelerations 
very  closely  and  were  not  included  in  any  of  the  inertial  loads  removal  from 
the  force/moment  measurements. 
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Testing  Ranges 

• Longitudinal  load  ranges 

(NF:  1000  lbs.;PM:  4000  in.-lbs.;  AF:  500  lbs.) 

• Accelerometer  range  (10  g’s) 

• Tunnel  velocity:  267  ft/sec 

• Model  sink  rates 

(0,  1,  2.33,  4.67,  7,  9.33,  1 1.667  ft/sec) 

• Gamma 

(0,  0.2, 0.5,  1,  1.5,2,  2.5  degrees) 

• Model  AOA  between  6 to  1 1 degrees 

• h/b  range 

(1.5  to  0.2) 


Here  is  the  range  of  variation  of  the  major  test  condition  parameters.  Note 
that  the  flight  path  angle  (Gamma)  corresponds  directly  to  model  sink  rate  for  a 
given  tunnel  velocity.  The  wind  tunnel  investigation  did  include  some  lower 
tunnel  velocity  test  points,  but  none  of  those  data  have  been  analyzed  yet  and 
are  not  included  in  this  presentation. 
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Data  Corrections 


• Static 

-weight  tare,  balance  woz 


• Dynamic 

-weight  tare,  balance  woz 

-model  acceleration  measurements  used  to  remove  inertial  loads  from 
balance  load  measurements 
-angle  of  attack  variation 


Typical  corrections  were  applied  to  the  static  data.  The  dynamic  data  included 
the  same  corrections  plus  those  necessary  to  remove  inertial  loads  and  any 
angle-of-attack  variation  during  the  plunge. 


2246 


Static  Ground  Effects 


• Art  Powell  will  present  the  next  section 
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Static  Ground  Effects  (SGE) 

• Flow  Angularity  (Model  10) 

• Lift  and  Moment  (Models  7,6) 

Objective:  Best  estimate  of  lift  and  pitching 
moment  characteristics  in  static  ground 
effect  for  comparison  with  dynamic  values. 


While  Lewis,  Sharon  Graves  (GWU  JIAFS  Graduate  Student),  and  Bob 
processed  the  dynamic  data.  Art  took  on  reducing  the  static  data  so  we  would 
have  a basis  of  comparison. 

Flow  angularity  was  measured  using  Model  10,  the  AR=7  NACA  0012 
elliptical  wing.  DACVINE,  a higher-order  panel  method,  was  used  to  estimate 
zero-lift  angle-of-attack,  since  this  model  had  an  underwing  fuselage  to  house 
the  balance. 

The  tunnel  with  DGE  strut  showed  a downflow  (of  approximately  0.15 
degrees)  at  the  tunnel  centerline,  which  decreased  to  near  zero  at  the  floor. 

The  flow  angularity  was  deemed  too  small  to  have  a significant  effect  on  the 
data  for  the  lower  AR  wings. 

Model  7,  the  Tu-144  planform,  and  Model  6,  the  TCA  planform,  were  the  only 
low-AR  Models  tested.  Severe  mechanical  problems  plagued  the  test 
throughout.  Data  quality  was  not  judged  good,  but  was  felt  to  be  sufficient  to 
determine  if  any  significant  DGE  effects  existed. 

Only  lift  and  moment  data  were  considered.  Since  the  models  had  no  camber, 
and  no  fuselage,  and  test  Reynolds  number  so  small,  there  was  not  sufficient 
motivation  to  study  the  drag  variation  with  ground  effect. 

Overall  agreement  of  static  data  with  other  sources  was  good. 
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Flow  Angularity  - Model  10 


Model  10  consisted  of  an  AR=7  elliptic  wing  with  symmetric  NACA0012 
sections  and  a fuselage  to  accommodate  the  balance.  Since  this  model  had  a 
steep  lift-curve  slope,  it  was  ideal  for  probing  tunnel  flow  angularity.  The 
presence  of  the  fuselage  made  the  model  slightly  asymmetric  with  respect  to 
lift,  so  DACVINE  was  used  to  determine  its  angle  of  zero  lift  at  various 
heights.  This,  along  with  angle  of  attack  sweeps  taken  at  various  heights,  was 
used  to  determine  the  flow  angularity.  The  figure  is  the  geometry  wireframe 
developed  for  the  DACVINE  analysis. 
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Static  Lift  Curves  for  Model  10 


Lift  data  were  taken  at  four  heights,  and  the  resulting  lift  curves  are  straight 
and  pass  close  to  the  origin.  The  abscissa  is  really  model  attitude  angle,  not 
angle  of  attack.  This  distinction,  which  is  usually  ignored,  is  necessary  in 
detecting  tunnel  flow  angularity. 


* 
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Zero-Lift  Offsets,  Model  10 


This  figure  is  an  expanded  view  near  the  origin  for  the  lift  curves  of  the  last 
figure.  As  can  be  seen,  the  curves  do  not  all  exactly  pass  through  the  origin. 
There  is  some  small  positive  angle  of  attack  indicated  at  zero  lift,  which 
changes  with  height  from  the  floor. 
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Tunnel  flow  angularity  is  calculated  using  the  zero-lift  angles  of  attack  from 
DACVINE,  and  the  zero-lift  attitude  angles  from  the  test.  The  blue  curve 
(solid  line)  is  the  resulting  flow  angularity  for  the  tunnel  with  this  cart  and 
model  in  place.  It  is  interesting  to  note  that  this  is  an  apparent  downflow,  the 
opposite  of  what  was  found  in  the  U&I  test  (LaRC442).  The  flow  angularity 
was  ignored  in  the  DGE  test  because  it  was  deemed  insignificant  for  the 
dynamic  measurements,  but  more  importantly  because  the  DGE  test  really  only 
required  that  a difference  be  seen  between  static  and  dynamic  data.  An  attempt 
is  underway  to  extract  “dynamic”  flow  angularity. 
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Model  7 Static  Characteristics 


Model  7 was  a flat  wing  of  Tu- 144  planform.  Its  leading  edge  was  of  small 
radius  from  centerline  to  tip.  It  featured  flaps  which  were  deflected  10  degrees 
for  this  test.  The  model  was  attached  to  the  sting  through  a balance  block  and 
balance  enclosed  in  a fairing  and  mounted  to  the  model’s  upper  surface.  A 9- 
degree  knuckle  was  used  to  provide  some  angle-of-attack  capability  near  the 
groundplane.  Model  size  is  indicated  above. 

The  model  was  run  through  angle-of-attack  sweeps  at  six  heights  above  the 
floor.  Floor  boundary  layer  suction  was  used. 
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Static  Lift  Curves,  Model  7 (Tu-144) 


6 7 8 9 10  11  12 

a (deg) 


Model  7 static  lift  characteristics  were  measured  at  a number  of  discrete 
heights,  which  for  comparison’s  sake  have  been  normalized  by  the  model  span, 
47.1”.  Each  data  point  represents  an  average  of  1200  samples,  which 
eliminates  sting  vibration,  a spurious  effect  present  in  most  of  the  data.  The 
data  are  shown  in  the  region  of  interest  to  the  DGE  test,  and  landing  maneuvers 
in  general.  The  h/b=.  125  data  is  limited  in  angle  of  attack  due  to  clearance 
problems  with  the  tunnel  floor. 
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Regressing  linearly  the  lift  data  in  the  7-1 1 degree  range  allows  the  lift  data  to 
be  generalized  and  interpolated,  and  although  it  was  not  strictly  necessary  in 
constructing  the  static  ground  effects  characteristics,  it  is  useful  for  data 
control.  The  lift  at  zero  angle  of  attack  asymptotes  at  higher  h/b,  which  is 
reassuring. 
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Regression  coefficients  were  quite  high  for  the  Model  7 lift  data.  Again,  the 
lift-curve  slope  asymptotes  at  high  h/b,  as  expected.  The  increase  in  lift  curve 
slope  at  low  h/b  values  is  substantial. 
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The  regressed  CL  @ a =0  and  dCL/da  vs.  h/b  curves  were  used  to  construct 
curves  of  CL  vs.  h/b  for  7,  9 and  1 1 degrees  angle  of  attack.  These  were 
normalized  by  the  free-air  CL  and  are  presented  here  as  normalized  ground 
effect  lift  increment.  The  data  compress  well  under  the  normalization.  The 
data  from  this  test  were  compared  with  data  from  the  Tu-144  flight  test,  in 
which  a series  of  landings  were  conducted  at  constant  angle  of  attack.  The 
agreement  is  very  good,  despite  the  fact  that  the  Tu-144  data  are  for  a dynamic 
maneuver  of  a full  configuration,  at  vastly  different  Reynolds  number,  and 
powered.  Also  note  the  large  (42%)  ground  effect  lift  increment  for  this 
planform. 
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Pitching  moment  data  about  the  50%-MAC  point  are  shown.  The  large 
positive  slope  indicates  that  the  aerodynamic  center  for  this  angle-of-attack 
range  is  forward  of  the  moment  reference.  The  data  do  show  a negative  (nose- 
down)  pitching  moment  increment  as  the  groundplane  is  approached,  which  is 
what  one  would  expect. 
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Re-referring  the  moment  data  to  42.5%-MAC  essentially  zeroes  the  moment 
slope,  indicating  this  to  be  the  approximate  free-air  aerodynamic  center  for  this 
angle-of-attack  range.  Thus  referred,  the  negative  moment  increment  due  to 
groundplane  proximity  is  larger  than  in  the  previous  figure.  The  normalized 
moment  increment  presented  in  the  next  slide  will  be  based  on  this  reference 
point. 
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This  figure  shows  the  normalized  moment  increment  about  the  aerodynamic 
center  for  Model  7.  Following  Curry,  the  ground-effect  moment  increment  is 
normalized  by  the  out-of-ground-effect  (OGE)  lift  coefficient.  The  data 
compresses  well  under  this  normalization.  One  Tu-144  point  is  shown.  This 
was  taken  from  an  eyeball  fairing  of  rather  noisy  flight  data,  an  indication  of 
the  noise  level  is  shown  by  the  error  bar.  Again,  the  agreement  appears  good, 
despite  likely  different  flap  settings. 


Reference: 

R.  E.  Curry,  “Dynamic  Ground  Effect  for  a Cranked  Arrow  Wing  Airplane,” 
NASA  TM-4799,  August  1997. 
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Model  6 Static  Characteristics 


T 


Model  6 was  the  TCA  planform.  It  had  no  camber  or  twist,  but  featured  10- 
degree  deflected  trailing  edge  flaps.  The  leading  edge  radius  was  as  large  as 
the  wing  thickness  allowed  inboard  of  the  leading-edge  sweep  break,  and  much 
smaller  outboard.  Like  Model  7 it  mounted  to  the  sting  via  a faired  balance 
block  and  balance  on  the  upper  surface.  The  9-degree  knuckle  was  used  with 
this  model  as  well. 

This  model  was  run  through  angle-of-attack  sweeps  at  five  heights  above  the 
floor.  Tunnel  boundary  layer  suction  was  used.  In  addition,  static  data  was 
available  at  high  h/b  from  certain  dynamic  runs  before  the  plunges  began. 
These  data  were  also  used  to  construct  an  estimate  of  the  aerodynamic 
characteristics  in  static  ground  effect. 
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Static  Lift  Curves,  Model  6 (TCA) 


« -TCA-1  (h/b=.174) 
-TCA-1  (h/b=. 987) 


The  Model  6 lift  curves  are  presented  in  the  angle  of  attack  range  of  interest  for 
DGE  testing.  Comparable  lift  curves  from  the  TCA-1  (LaRC  14’x22’  Test 
449)  are  shown  also.  The  h/b=.987  (TCA-1)  data  (Run  572)  should  be 
considered  as  free  air  data,  since  it  includes  wall  corrections.  The  h/b=.174 
(gear  height)  data  (Run  350)  was  processed  as  ground  effects  data  and  was  not 
corrected  for  the  floor  presence.  Both  of  the  TCA-1  curves  are  for  30  deg  TE 
flaps,  as  opposed  to  the  10-degree  TE  flaps  of  the  DGE  model.  The  TCA-1 
curves  are  for  a full  configuration,  but  without  horizontal  tail,  while  the  DGE 
test  data  are  for  a flat  wing  planform  only.  The  Reynolds  numbers  are  also 
different:  7.8  million  for  the  TCA-1  data  versus  only  4.7  million  for  the  DGE 
test  data.  The  h/b=.174  TCA-1  data  was  limited  to  below  about  8.6  degrees  by 
tailstrike. 

The  bulk  of  the  DGE  test  data  came  from  constant-height  angle-of-attack 
sweeps.  For  these  points,  roughly  1200  data  samples  are  averaged  for  each 
data  point.  Data  at  h/b=1.87  came  from  DGE  runs  prior  to  plunge  start.  Lift- 
curve  data  are  available  because  plunges  were  taken  at  7, 9 and  1 1 degrees 
angle  of  attack  for  Model  6.  Typically,  over  100  data  samples  were  available 
for  time-averaging  prior  to  plunge  start.  This  largely  eliminates  the  effects  of 
sting  oscillations,  if  any. 

A few  individual  data  points  at  low  h/b  are  shown  near  7,  9,  and  1 1 degrees 
angle  of  attack.  These  points,  while  unique  in  h/b,  have  turned  out  to  be 
extremely  useful,  as  will  be  seen. 
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Each  of  the  DGE  test  static  lift  curves  was  linearly  regressed  between  7 and  1 1 
degrees.  Regression  coefficients  “r”  typically  were  quite  good,  with  at  least 
three  “nines”  past  the  decimal.  The  lift  curves  were  represented  as  a zero-a  CL 
and  lift  curve  slopes.  The  zero-a  CL  data  are  presented  here  as  a function  of  h / 
b.  These  data  appear  well-behaved  and  plausible. 
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This  plot  shows  the  regressed  lift-curve  slopes  of  the  Model  6 static  ground 
effects  data.  One  expects  this  curve  to  asymptote  to  some  value  as  h/b 
approaches  2,  which  it  does  if  the  data  point  at  h/b=1.5  is  ignored.  In  the 
analysis,  lift  data  at  this  h/b  is  excluded  in  favor  of  the  h/b=1.87  data  taken 
before  plunge  start  on  the  dynamic  runs. 
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With  the  regressed  data  of  the  last  two  figures  it  is  possible  to  put  together 
plots  of  ground  effect  CL  increment  versus  h/b  at  constant  angle  of  attack. 

This  is  a plot  of  that  data,  normalized  by  free-air  CL.  Since  the  lowest  h/b  for 
which  a lift  curve  could  be  constructed  was  0.250,  extension  to  gear  height 
might  seem  difficult,  except  that  the  extra,  unique  data  points  at  lower  h/b  were 
very  close  to  the  angles-of-attack  for  which  the  plot  was  made.  A short 
extrapolation,  using  extrapolated  lift-curve  slopes  makes  the  low  h/b  end  of  the 
curve  accessible.  The  TCA-1  data  are  shown  here  for  comparison,  and  show 
reasonable  agreement  with  the  DGE  static  data.  The  9-deg  TCA  data  required 
a short  extrapolation,  since  the  model  was  tailstrike  limited  to  8.6  degrees  at  h/ 
b=.174. 

Note  that  the  static  lift  increment  in  ground  effect  for  this  planform  is 
significantly  lower  (-25-30%  compared  to  42%)  than  it  is  for  the  Tu-144. 
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This  is  a plot  of  the  pitching  moment  coefficient  taken  about  the  50%  MAC  for 
Model  6.  At  this  moment  reference  location,  decreasing  h/b  causes  a positive 
Cm  increment.  This  and  the  positive  slope  suggests  that  the  moment-reference 
center  is  well  behind  the  aerodynamic  center  for  this  angle-of-attack  range. 
Also,  note  that  the  h/b=1.87  data  lie  well  below  the  other  curves,  which  is 
counterintuitive.  One  would  expect  the  curves  to  be  close  together  at  high  h/b, 
near  the  edge  of  the  ground  effect  regime. 
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The  moment  data  was  re -referenced  to  35.5%  MAC,  which  reduced  dCm/da  to 
essentially  zero  for  h/b=  1.5.  Taken  about  this  reference  point,  most  of  the  data 
exhibit  a negative  moment  increment  as  the  groundplane  is  approached.  The 
data  at  h/b=1.87,  taken  from  the  dynamic  runs  prior  to  plunge  start,  do  not 
follow  this  trend,  which  is  counterintuitive  and  at  odds  with  the  other  data. 
These  data  are  therefore  not  used  for  OGE  reference  in  normalizing  the 
moment  increment  data. 


2267 


HSCTHiffh  Lift  Aerodynamics 


This  figure  shows  the  normalized  moment  increment  due  to  static  ground 
effect  as  measured  for  Model  6.  The  data  show  the  expected  trend  of  negative 
(nose-down)  pitching  moment  increment  as  the  groundplane  is  approached. 
The  data  does  not  collapse  as  well  under  the  normalization  as  did  those  for 
Model  7.  Limited  data  from  TCA-1  appears  to  agree  with  the  trend  shown  by 
these  data. 
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Static  Ground  Effect 
Conclusions 

• Flow  angularity  check  shows  small  downwash, 
which  was  ignored  in  the  DGE  analysis. 

• Model  7 data  agrees  with  Tu-144  flight  test. 

• Model  7 lift  and  moment  data  compressed  well 
under  normalization. 

• Model  6 data  compares  well  with  TCA- 1 . 

• Model  6 moment  data  did  not  compress  well 
under  normalization. 


The  static  data  serve  as  a basis  for  comparison  for  the  dynamic  data.  An 
examination  of  the  tunnel  flow  angularity  found,  for  the  DGE  strut  and  Model 
10,  a slight  downwash  of  about  .15  degrees  near  the  tunnel  centerline,  which 
decreased  as  the  tunnel  floor  was  approached.  Since  the  DGE  test  was 
principally  a comparison  of  static  and  dynamic  effects,  and  based  on  the 
assumption  that  flow  angularity  would  be  the  same  for  dynamic  conditions  as 
for  static  conditions,  the  flow  angularity  was  ignored  for  the  purposes  of  this 
study. 

The  Model  7 static  lift  and  moment  data  was  quite  consistent,  agreed  well  with 
data  from  Tu-144  flight  test,  and  collapsed  well  under  normalization. 

Model  6 static  data  was  less  consistent,  but  lift  and  moment  increment  data 
agreed  reasonably  well  with  TCA-1  wind-tunnel  data.  The  lift  data  collapsed 
reasonably  well  under  normalization,  but  the  moment  data  did  not. 
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Dynamic  Ground  Effects 


• Lewis  Owens  will  present  the  next  section 
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Dynamic  Ground  Effects  (DGE)  Data 
Challenges 

• Different  sampling  rates 

-model  height  & pitch  position  (50  to  60  hertz) 

-balance  and  accelerometers  (150  hertz) 

• Accelerations  not  completely  zeroed 

• Model/support  system  “ringing” 

-during  constant  velocity  segment 

-most  noticeable  at  high  sink  rates  (>4  ft/sec) 


Some  aspects  of  the  dynamic  data  acquisition/reduction  presented  real  challenges.  The  data 
system  was  set  up  to  acquire  multiple  channels  of  data  that  were  each  scanned  at  150  hertz  over 
an  8 second  sampling  period.  Some  of  the  more  significant  problems  with  the  dynamic  data 
included  issues  associated  with  sampling  rates,  zeroing  initial  accelerations  and  flexibility  of 
the  model  support  system.  The  position  measurements  were  not  sampled  at  as  high  a rate  as  the 
balance  and  acceleration  measurements.  This  problem  is  obvious  when  either  balance  or 
accelerometer  data  is  plotted  versus  height.  Multiple  data  points  are  acquired  at  a given  height 
measurement.  The  impact  of  this  on  the  DGE  plots  is  not  considered  significant  for  the  height 
measurements,  which  results  in  small  shifts  in  the  data  on  the  order  of  0.01  to  0.07  in  h/b 
depending  on  the  sink  rate  of  the  run.  The  impact  of  reduced  sampling  rate  on  the  pitch 
measurement  is  not  as  clear  cut  and  may  be  contributing  to  problems  with  cleaning  up  the  rest 
of  the  spread  in  the  force/moment  data  already  corrected  for  inertial  loads.  No  corrections  for 
these  sampling  rate  problems  have  been  made  to  any  of  the  dynamic  data  presented.  Another 
challenge  included  the  initial  zeroing  of  the  model  acceleration  levels  at  the  beginning  of  each 
dynamic  plunge.  We  attempted  to  handle  this  with  a combination  of  wind-off  zero  and 
theoretical  corrections  for  changes  in  the  gravity  component  with  pitch,  but  this  did  not  work. 
Sharon  Graves  and  I ended  up  taking  an  average  level  at  the  beginning  of  each  run  to  reference 
the  accelerations.  This  appeared  to  work  very  well  and  allowed  us  to  use  the  integrated 
accelerations  to  calculate  the  model’s  sink  rate.  Finally,  model  support  system  flexibility 
tended  to  shorten  the  constant  velocity  segment  of  the  dynamic  plunge  trajectories  and  were 
associated  with  the  larger  spread  in  the  data  for  higher  sink  rate  runs.  In  the  future,  we  plan  for 
changes  in  the  cart  control  system  to  help  alleviate  this  situation. 
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Sample  of  Inertial  Load  Removal 
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Here  is  a representative  example  of  the  effect  of  inertial  load  removal  from  the 
normal  force  balance  measurements  during  a lower  sink  rate  run.  A run  begins 
with  the  model  resting  at  a position  near  the  tunnel  centerline,  in  this  case  at  an 
h/b  value  of  about  1.90.  As  the  dynamic  plunge  starts,  large  excursions  in  the 
normal  load  is  evidence  of  the  rapid  acceleration  to  the  target  sink  rate.  At  an 
h/b  level  of  about  1.30,  the  constant  velocity  segment  of  the  trajectory  is 
reached  and  the  dynamic  load  excursions  have  diminished.  The  steady 
increase  in  the  normal  force  level  as  h/b  decreases  is  associated  with  the 
ground  effect.  At  an  h/b  of  about  0.3,  the  model  begins  decelerating  and  the 
raw  normal  force  level  is  consistently  lower  than  the  corrected  data.  The  dash 
curve  represents  the  normal  force  data  after  corrections  for  the  primary  inertial 
loads.  The  amount  of  inertial  load  clean  up  is  most  evident  in  comparing  the 
raw  and  primary  curves  near  an  h/b  of  1 .4  and  below  an  h/b  of  0.3.  Between 
these  two  h/b  values,  the  primary  inertial  load  corrections  routine  does  not 
have  a significant  impact  on  the  raw  data. 
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Sample  of  Inertial  Load  Removal 
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Here  is  a representative  example  of  the  effect  of  inertial  load  removal  from  the 
normal  force  balance  measurements  during  a higher  sink  rate  ran.  For  the 
higher  sink  rate  runs,  the  effect  of  removing  the  primary  inertial  loads  is  more 
evident  throughout  the  plunge  trajectory.  Notice  that  the  model  deceleration 
starts  earlier  (h/b  level  of  about  0.4)  as  compared  with  the  previous  plot  for  the 
lower  sink  rate  run.  Also  note  that  the  support  system  “ringing”  is  more 
prevalent  throughout  the  plunge. 
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Trajectory  Parameter:  Model  Gamma 
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Here  is  a representative  plot  of  the  model  glide  path  angle  that  resulted  for  the 
different  dynamic  plunge  trajectories.  The  model  sink  rates  were  chosen  to 
provide  selected  constant  gamma  for  the  given  tunnel  velocity  (267  ft/sec).  A 
significant  feature  of  these  trajectories  is  the  long  run  of  constant  gamma  for 
the  lower  sink  rate  runs.  Notice  that  this  constant  gamma  segment  becomes 
shorter  with  increasing  sink  rate  because  the  model  support  ringing  persists  to 
lower  values  of  h/b.  Also  note  that  the  model  deceleration  starts  earlier  (higher 
h/b  values)  as  the  sink  rate  increases,  which  further  restricts  or  eliminates  the 
constant  gamma  segment.  Another  noticeable  feature  in  this  data  is  associated 
with  the  reduced  height  sampling  rate.  For  the  sink  rates  greater  than  1 ft/sec, 
the  data  symbols  tend  to  cluster  in  groups  of  three  showing  that  three 
accelerometer  measurements  were  acquired  before  the  height  measurement 
was  updated. 
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Trajectory  Parameter:  Model  Alpha 
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Here  is  a representative  plot  of  the  model  angle  of  attack  that  resulted  for  the 
different  dynamic  plunge  trajectories.  As  with  the  glide  slope  angle,  the 
constant  model  alpha  segment  was  longer  for  the  lower  sink  rate  runs.  For  the 
higher  sink  rate  runs,  the  model  alpha  varied  more  significantly.  Since  the 
main  point  of  this  investigation  was  to  try  and  isolate  the  ground  effect  (that  is, 
look  at  lift  variation  with  constant  alpha),  we  decided  to  correct  the  lift  data  to 
a constant  alpha  to  remove  this  effect.  This  correction  involved  compensating 
for  the  alpha  variation  by  using  the  out-of-ground  effect  (OGE)  lift  curve 
slope.  All  of  the  DGE  data  that  will  be  presented  has  been  corrected  to  a 
constant  alpha  to  make  it  comparable  to  the  SGE  data. 
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Model  7 DGE  Results:  Lift 
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Here  is  a plot  of  the  lift  increment  due  to  ground  effect  for  the  dynamic  data 
obtained  for  Model  7.  The  cluster  of  data  points  at  an  h/b  of  1.9  gives  an  idea 
of  the  variation  in  the  lift  levels  while  sitting  statically  before  the  start  of  the 
dynamic  plunge.  The  data  during  the  dynamic  plunge  tends  to  cluster  in  a 
solid  trend  indicating  the  ground  effect  levels  on  lift.  The  data  points  falling 
outside  this  trend  are  associated  with  the  higher  sink  rate  data.  Keep  in  mind 
that  up  to  h/b  levels  of  0.4,  the  model  is  decelerating  for  the  higher  sink  rate 
runs  so  that  the  sink  rate  is  not  constant.  From  this  DGE  lift  data,  there  is  no 
indication  of  a significant  change  in  the  ground  effect  trend  for  varying  sink 
rates. 
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Here  is  a plot  of  the  pitching  moment  increment  due  to  ground  effect  for  the 
dynamic  data  obtained  for  Model  7.  The  cluster  of  data  points  at  an  h/b  of  1 .9 
gives  an  idea  of  the  variation  in  the  pitching  moment  levels  while  sitting 
statically  before  the  start  of  the  dynamic  plunge.  The  data  during  the  dynamic 
plunge  tends  to  cluster  in  a solid  trend  indicating  the  ground  effect  levels  on 
pitching  moment.  The  ground  effect  trend  on  normalized  pitching  moment 
magnitude  is  about  2 percent.  If  this  data  were  referenced  to  a moment  center 
closer  to  the  aerodynamic  center  (42%  mac),  then  the  magnitude  of  this  effect 
would  increase  to  about  5 percent. 
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Model  7 DGE/SGE  Results:  Lift 
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Since  there  did  not  appear  to  be  any  significant  difference  in  the  ground  effect 
trend  with  different  model  sink  rates,  one  of  the  lower  sink  rate  DGE  data  runs 
was  selected  for  comparison  with  the  SGE  data.  This  run  was  chosen  because 
it  was  the  highest  sink  rate  run  available  in  which  the  model  deceleration  did 
not  begin  until  it  was  below  an  h/b  of  0.3.  This  permits  a comparison  of  DGE 
lift  increase  levels  with  those  of  the  SGE  in  the  more  sensitive  region  of  the 
ground  effect  trend.  This  comparison  does  not  show  that  the  DGE  lift 
increase  is  significantly  different  from  that  of  the  SGE.  (Recall  that  we  are 
looking  for  differences  indicating  that  the  DGE  lift  increment  is  about  50 
percent  of  that  of  the  SGE  at  values  of  h/b  of  about  0.3.) 
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Model  6 DGE  Results:  Lift 
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Here  is  a plot  of  the  lift  increment  due  to  ground  effect  for  the  dynamic  data 
obtained  for  Model  6.  Again,  the  cluster  of  data  points  at  an  h/b  of  1.9  gives 
an  idea  of  the  variation  in  the  lift  levels  while  sitting  statically  before  the  start 
of  the  dynamic  plunge.  The  data  during  the  dynamic  plunge  also  tends  to 
cluster  in  a solid  trend  indicating  the  ground  effect  levels  on  lift.  The  data 
points  falling  outside  this  trend  are  associated  with  the  higher  sink  rate  data. 
Keep  in  mind  that  up  to  h/b  levels  of  0.4,  the  model  is  decelerating  for  the 
higher  sink  rate  runs  so  that  the  sink  rate  is  not  constant.  From  this  DGE  lift 
data,  there  is  no  significant  change  in  the  ground  effect  trend  for  varying  sink 
rates.  Note  that  the  ground  effect  lift  increase  for  this  wing  planform  is  below 
15  percent  at  an  h/b  of  0.3  as  compared  to  about  20  percent  for  Model  7.  Also, 
the  lift  increment  data  for  Model  6 appears  to  have  a larger  variation  band  than 
that  for  Model  7.  Model  6 had  a blunt  inboard  LE  radius  and  Model  7 had  a 
“sharp”  LE.  This  geometry  difference  may  contribute  to  the  larger  variation 
in  the  data  for  Model  6. 
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Here  is  a plot  of  the  pitching  moment  increment  due  to  ground  effect  for  the 
dynamic  data  obtained  for  Model  6.  The  cluster  of  data  points  at  an  h/b  of  1 .9 
gives  an  idea  of  the  variation  in  the  pitching  moment  levels  while  sitting 
statically  before  the  start  of  the  dynamic  plunge.  The  data  during  the  dynamic 
plunge  was  more  scattered  than  that  for  Model  7.  The  ground  effect  trend  on 
the  pitching  moment  increment  is  not  clear. 
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Model  6 DGE/SGE  Results:  Lift 


Since  there  did  not  appear  to  be  any  significant  difference  in  the  ground  effect 
trend  with  different  model  sink  rates,  one  of  the  lower  sink  rate  DGE  data  runs 
was  selected  for  comparison  with  the  SGE  data.  This  run  was  chosen  to  be 
consistent  with  that  chosen  for  Model  7.  This  comparison  also  does  not  show 
that  the  DGE  lift  increase  is  significantly  different  from  that  of  the  SGE. 
(Recall  that  we  are  looking  for  differences  indicating  that  the  DGE  lift 
increment  is  about  50  percent  of  that  of  the  SGE  at  values  of  h/b  of  about  0.3.) 
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Why  are  SGE/DGE  similar? 


• Previous  data  indicated  DGE/SGE  lift  increment 
ratio  of  about  50%  for  aspect  ratios  less  than  2.0 

• LE  sweep  or  wing  sweep-related  factor  may  be  the 
more  important  controlling  parameter  (especially 
with  discontinuous  LE) 

• Other  factors  may  include:  AOA  region,  RN,  flap 
configuration 


The  14x22  ft  DGE  test  did  not  show  significant  differences  between  SGE  and 
DGE  lift  increment  data.  The  previous  KU  data  suggests  that  the  models 
tested  were  in  the  aspect  ratio  range  where  this  difference  should  be 
significant.  A review  of  the  KU  database  was  performed  to  check  for 
consistency  and  what  factors  may  explain  the  difference  in  the  findings.  One 
difference  noted  was  that  the  models  tested  in  the  14x22  ft  study  had  breaks  in 
the  wing  planform  LE  while  the  KU  data  was  based  on  planforms  with  no  LE 
breaks.  For  continuous  LE  planforms,  there  is  a direct  correlation  between  the 
LE  and  the  aspect  ratio.  For  planforms  with  a LE  break,  this  relationship  is  not 
as  direct.  Other  differences  were  also  noticed  and  considered  in  the  review  of 
the  ground  effects  database.  However,  the  LE  sweep  factor  seemed  to  be 
consistent  and  deserved  further  examination. 


University  of  Kansas  data  reference: 

R.  C.  Chang  and  V.  U.  Muirhead,  “Effect  of  Sink  Rate  on  Ground  Effect  of 
Low-Aspect-Ratio  Wings,”  Journal  of  Aircraft,  Vol.  24,  No.  3,  March  1987, 
pp.  176-180. 
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Here  is  another  look  at  the  University  of  Kansas  wind-tunnel  data  plotted 
versus  LE  sweep  (outboard  LE  for  cranked  wing  planforms)  instead  of  aspect 
ratio.  Comparable  SGE  and  DGE  data  points  were  added  to  show  where  they 
fall  relative  to  the  existing  data.  Note  that  both  the  TCA  and  the  Tu-144  model 
data  was  placed  considering  the  sweep  of  the  outboard  LE.  If  this  is  a proper 
way  of  looking  at  the  differences  between  SGE/DGE  data,  then  this  may 
explain  why  the  models  tested  in  the  14’x22’  study  did  not  show  any 
significant  difference. 
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Computational  Checks 

• Unsteady  3-D  seems  to  be  predicting  the  similar 
DGE  results  as  experiment  for  a 60  degree  delta 
wing  and  the  TCA  wing 

• Recommend  using  computational  methods  to 
investigate  other  factors  (LE  sweep,  AGA 
region,...) 


Two  different  analytical  approaches,  each  developed  by  Winfried  Feifel 
(Boeing)  and  Bill  Dwyer  (Northrup  Grumman),  were  used  to  predict  dynamic 
ground  effects.  Both  the  analytical  and  the  experimental  DGE  results  are 
consistent  and  show  no  significant  difference  between  SGE  and  DGE  for  the 
TCA  wing  planform.  It  is  interesting  to  note  that  Bill  Dwyer  showed  some 
potentially  significant  DGE/SGE  differences  for  the  XB-70  wing,  but 
unfortunately  this  work  was  not  concluded  due  to  funding  constraints. 


These  analytical  approaches  provide  an  opportunity  to  explore  other  factors 
associated  with  low  aspect  ratio  wings  to  gain  a better  understanding  of  when 
there  may  be  significant  differences  between  SGE  and  DGE  testing  techniques. 
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Dynamic  Ground  Effects:  Flight  Data 

• Bob  Curry  will  present  the  next  section 
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Flight  Data 

• Flight  data  obtained  for  two  similar  vehicles 

- F- 1 6 XL  experiment  complete 

- TU- 144,  preliminary  data  available 

• Data  obtained  Test  Techniques 

- Instrumentation  included  both  on-board  sensors  and  differential  GPS 

- Data  was  obtained  during  constant-angle-of-attack,  constant-thrust 
approaches  to  the  runway 

• Data  Analysis 

- Data  corrected  for  variations  in  alpha  and  surface  deflections 

- Flight  results  are  shown  as  untrimmed,  incremental  changes  in  lift, 
and  pitching  moment  coefficient  relative  to  the  ‘out-of-ground- 
effect’  lift  coefficient 


A flight  experiment  was  conducted  to  provide  additional  information  regarding 
ground  effect  characteristics  for  slender-wing,  high-speed  configurations. 
Flight  data  has  been  obtained  for  the  F-16XL  and  the  Tupolev  TU-144 
supersonic  aircraft.  Data  from  the  TU-144  flight  experiment  will  also  be  used 
to  validate  results  from  the  DGE  wind  tunnel  test  technique. 


A thorough  discussion  of  the  flight  test  techniques  and  data  analysis  process 
for  the  F-16XL  is  provided  in  NASA  TM  4799.  Similar  methods  are  being 
used  for  the  TU-144  flight  experiment. 


The  data  were  obtained  during  approaches  to  the  runway.  Before  each 
maneuver,  the  pilot  began  a stabilized  descent  at  a pre-determined  glide  slope 
and  angle  of  attack.  The  pilot  attempted  to  hold  the  power  and  angle-of-attack 
constant  for  the  remainder  of  the  maneuver  as  the  airplane  descended  into 
ground  effect.  Perturbations  from  the  initial  flight  conditions  were  attributed 
to  ground  effect.  Adjustments  to  angle-of-attack  and  elevon  position  which 
occurred  during  the  maneuvers  were  accounted  for  in  the  extraction  of  ground 
effect  increments. 


Data  were  obtained  for  a variety  of  gross  weights,  flight  path  angles,  flap  and 
canard  positions. 
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This  figure  shows  three  configurations  for  which  static  and  dynamic  ground 
effect  data  have  been  obtained.  All  of  the  wings  have  similar  aspect  ratio  and 
similar  inboard  and  outboard  leading  edge  sweep. 


The  wind  tunnel  model  is  an  uncambered,  planform  of  the  TU-144  wing.  The 
trailing  edge  elevons  were  set  at  10  degrees  (trailing  edge  down)  for  all  tests. 


Flight  data  for  the  TU-144  was  obtained  using  a research  vehicle  derived  from 
an  early  production  supersonic  transport.  The  configuration  has  four  turbojet 
engines  mounted  below  the  wing  similar  to  ‘HSCT’  configurations.  The 
airplane  is  normally  landed  with  canard  deployed,  however,  flight  test  data  was 
obtained  with  the  canard  both  deployed  and  retracted.  Elevon  positions  were 
were  about  9 deg  (trailing  edge  down)  during  final  approach  for  the  majority  of 
flight  maneuvers.  Slightly  negative  elevon  positions  (trailing  edge  up)  were 
required  during  maneuvers  with  the  canard  retracted  . 


The  F-16XL  airplane  is  a high-performance,  single-seat  airplane  with  a 
cranked-arrow-wing  designed  for  supersonic  cruise  flight.  The  single  turbofan 
engine  is  located  on  the  centerline.  Because  of  the  relatively  large  gear  height 
relative  to  wing  span,  ground  effect  data  at  very  low  ratios  of  h/b  could  not  be 
obtained  with  the  F- 1 6XL. 


2287 


During  flight  testing,  the  flight  path  angle  and  vertical  speed  is  generally 
constant  as  the  airplane  descends  through  altitudes  above  ground  effect. 

Initial  flight  path  angles  for  the  flight  maneuvers  varied  between  about  -2  and 
-3  deg.  These  sink  rates  correspond  to  vertical  velocities  of  between  10  and 
15  ft/sec.  As  the  airplane  passes  through  a height  of  about  one  span,  the 
flight  path  angle  and  sink  rate  begins  to  decrease  due  to  ground  effect,  until 
both  parameters  are  nearly  zero  at  touchdown 


The  variation  of  flight  path  angle  with  altitude  during  the  flight  test  maneuvers 
is  similar  to  the  landing  flare  conducted  during  normal  landings  of  the  TU-144. 

The  DGE  wind  tunnel  capability  allowed  data  to  be  collected  at  a constant 
flight  path  angle  in  the  presence  of  ground  effect.  Similar  results  have  not 
been  obtained  in  flight  testing.  The  DGE  wind  tunnel  flight  path  angles  are 
generally  lower  than  the  sink  rates  experienced  in  flight. 
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During  the  flight  test  maneuvers,  the  pilot  attempts  to  hold  angle-of-attack 
constant  during  the  descent  through  ground  effect.  Atmospheric  disturbances 
and  the  need  to  insure  satisfactory  conditions  at  touchdown  cause  some 
variations  to  occur.  For  most  of  the  TU-144  flight  maneuvers  analyzed  so  far, 
there  has  been  a decrease  in  alpha  just  prior  to  touchdown.  Angle-of-attack  is 
generally  increased  slightly  during  conventional  landings.  During  the  flight 
test  maneuvers,  the  angle  of  attack  was  about  9 degrees. 

Angle-of-attack  was  constant  during  a significant  range  of  the  DGE  wind 
tunnel  trajectories.  As  the  model  decelerated  at  the  end  of  a run,  the  angle-of- 
attack  decreased.  The  region  of  constant  angle-of-attack  was  therefore  smaller 
for  the  higher  sink  rate  tests. 

Analysis  of  both  the  wind  tunnel  and  flight  data  sets  account  for  any  variations 
of  angle-of-attack.  Static  measurements  of  lift  and  pitching  moment 
derivatives  with  respect  to  angle-of-attack  were  used  for  these  corrections. 

Similar  corrections  are  made  to  the  flight  data  for  any  variation  in  elevon 
position  which  occurs  during  a maneuver. 

As  a result,  the  final  wind  tunnel  and  flight  data  sets  represent  untrimmed  force 
and  moment  increments  due  to  ground  effect. 
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The  incremental  changes  in  lift  coefficient  due  to  ground  effect  from  several 
flight  maneuvers  are  repeatable  within  the  noise  level  of  the  measurements. 
The  primary  sources  of  scatter  in  the  data  are  felt  to  be  noise  in  the 
accelerometer  measurements  and  atmospheric  disturbances. 


Data  for  both  the  TU-144  and  F-16XL  compare  favorably,  despite  differences 
in  wing  shape  and  engine  configuration. 
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Preliminary  comparisons  of  flight  measured  lift  coefficient  increments  from 
flight  and  DGE  wind  tunnel  testing  show  excellent  agreement.  The  flight  data 
shown  was  obtained  at  an  initial  sink  rate  of  approximately  15  ft/sec  with  the 
canard  deployed.  The  DGE  data  was  obtained  at  a sink  rate  of  about  2.33  ft/ 
sec,  wing  planform  only. 

Additional  correlations  will  be  possible  as  more  TU-144  flight  data  becomes 
available. 
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HSCT High  Lift  Aerodynamics 


Comparison  of  Wind  Tunnel/Flight  Pitching 


0.05 

0.04 

o 

CB 

E 0.03 


ge 

9 

VI 

<§> 

<y 

oc 

0 

JT 

u 

1 

u 


o> 

Q 


0.02 

0.01 

0.00 

-0.01 

-0.02 

-0.03 

-0.04 

-0.05 


c 

: ° 

m 

O Tu-144  Plan  form:  -2.33  Wstc 
— •—  Tu-144:  Flight  4 

OO 

■ 

n 

■ 

IS 

■ 

Q 

■a 

iSi 

V -Ms?! 

iS 

s yjj 

0 

u 

SESI 

3 

r 

m 

■a 

ms 

m 

11 

c 

m 

s 

■i 

0.00  0.20  0.40  0.60  0.80 


1.00 

h/b 


1.20  1.40  1.60  1.80  2.00 


Both  flight  and  wind-tunnel  data  indicated  small  levels  of  nose-down  pitching 
moment  in  the  presence  of  ground  effect  when  referenced  to  50%  MAC. 


The  center  of  gravity  for  the  aircraft  is  typically  about  40%  MAC  and  therefore 
the  negative  pitching  moment  due  to  ground  effect  is  more  significant  for  the 
flight  vehicle. 


Any  differences  between  the  wind  tunnel  and  preliminary  flight  data  are  within 
the  noise  level  of  the  flight  measurements. 
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HSCT  High  Lift  Aerodynamics 

Status  of  Flight  Testing 

• TU-144  Flight  Testing  Still  in  Progress 

• 17  flights  completed  to  date,  9 ground  effect 
maneuvers  attempted 

• Post-flight  computations  of  mass  properties  and 
thrust  estimates  for  most  maneuvers  still  in 
progress 

• Flight  data  set  to  be  summarized  in  TU-144  final 
report 


The  TU-144  flight  experiment  is  still  in  progress.  A total  of  18  flights  are 
expected  to  be  completed  by  the  end  of  February. 


The  flight  data  set  will  be  limited  to  about  8 to  10  good  maneuvers.  Although 
the  airplane  conducted  fifteen  flights,  ground  effect  maneuvers  were  only 
attempted  on  certain  flights  due  to  test  planning  constraints,  and  it  was  also  not 
possible  to  obtain  more  than  one  test  point  on  a given  flight. 


Additional  smoothing  and  trajectory  reconstruction  methods  may  be  used  to 
enhance  the  quality  of  the  flight  data.  Final  data  analysis  will  not  be  possible 
until  post-flight  computations  of  mass  properties  and  thrust  have  been  provided 
by  Tupolev.  A summary  of  the  flight  experiment  and  electronic  files 
containing  flight  data  (in  engineering  units)  will  be  available  to  the  HSR 
project  at  the  conclusion  of  the  flight  program. 
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Composite  Data  Set  as  a Function  of 


Shown  in  this  figure  is  the  composite  ground  effect  data  set  from  several  tests 
plotted  as  a function  of  aspect  ratio.  As  previously  noted,  there  are  two  general 
trends. 


One  trend  line  consists  of  all  the  dynamic  data  as  well  as  some  static  data. 
Another  trend  line  consists  of  several  static  test  points,  with  values 
approximately  twice  as  great  as  the  first  trend  line.  The  Tu-144  and  TCA 
dynamic  and  static  wind  tunnel  data,  as  well  as  the  Tu-144  flight  data  all  fall 
along  the  initial  trend  line.  The  dynamic  flight  data  for  the  F-16XL  also  fits 
this  trend,  however,  the  corresponding  F-16XL  static  wind  tunnel  data  is 
significantly  larger  than  the  flight  data. 
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When  the  same  data  is  plotted  as  a function  of  outboard  leading  edge  sweep 
angle,  both  the  Tu-144  and  TCA  data  from  the  current  study  fit  the  trends 
fairly  well.  The  F-16XL;  however,  appears  to  depart  from  the  trend.  It  should 
be  noted  however,  that  the  leading  edge  sweep  breakpoint  is  relatively  farther 
outboard  for  the  F-16XL  than  for  either  the  Tu-144  or  the  TCA  model. 
Therefore  it  could  be  argued  that  the  F-16XL  data  should  be  plotted  against  its 
inboard  leading  edge  sweep  angle  (70  deg),  which  would  result  in  a better  fit  to 
the  trend.  These  observations  suggest  that  ground  effect,  and  the  sensitivity  to 
dynamic  ground  effect,  may  be  indicated  by  a parameter  which  involves  a 
weighted  value  of  leading  edge  sweep. 


The  current  data  base  is  awkward  to  interpret  since  it  includes  an  inconsistent 
variety  of  configurations  and  test  facilities.  In  many  cases,  dynamic  data  for  a 
configuration  was  obtained  from  flight  or  a dynamic  facility  and  then 
compared  against  static  data  from  a different  wind  tunnel  facility.  Only  the 
University  of  Kansas  low  speed  wind  tunnel  and  the  NASA  Langley  14x22  ft 
DGE  cart  have  the  capability  to  provide  parametric  data  under  both  static  and 
dynamic  conditions.  So  far,  the  Kansas  data  have  distinctly  indicated  strong 
influences  while  the  14x22  ft  data  have  indicated  negligible  differences  due  to 
dynamics.  Further  parametric  testing  with  the  14x22  DGE  system,  especially 
with  higher  wing  sweep  angles,  may  help  isolate  the  controlling  parameters  for 
dynamic  ground  effects. 
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HSCT  High  Lift  Aerodynamics 


Conclusions 

• For  the  low  aspect  ratio  wing  planforms  tested 

- no  significant  DGE/SGE  differences 

- other  parameters  (besides  AR)  may  be 

significant  controlling  factors  in  this  difference  (LE  sweep) 


• Wind  tunnel  and  flight  ground  effects  increment 
data  for  the  Tu-144  compared  well 


For  the  HSCT  wing  planforms  tested  in  the  14x  22  ft  DGE  test,  no  significant 
differences  were  found  between  DGE  and  SGE  test  techniques.  From  previous 
ground  effects  data,  the  aspect  ratios  of  the  model  wing  planforms  tested  were 
such  that  differences  in  DGE  and  SGE  data  were  expected.  Closer 
examination  of  all  the  data  suggested  that  other  factors  (in  addition  to  AR)  may 
need  to  be  controlled  to  better  understand  this  difference.  Comparisons  of  the 
ground  effects  increment  data  from  the  14x22  ft  DGE  test  and  the  flight  test  for 
the  Tu-144  were  good.  These  ground  effects  increments  compared  well  even 
with  a very  basic  model  that  represented  only  the  wing  planform  of  the  Tu-144 
aircraft. 
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HSCT  High  Uft  Aerodynamics 

Recommendations 

• DGE  cart  needs  to  be  reworked  before  it  can  be  a 
routinely  used  in  ground  effects  testing 

• Follow-on  DGE  test 

- Check-out  reworked  DGE  cart 

- Confirm  current  HSR  DGE  findings  and  provide  an  opportunity  to 
validate  DGE  test  technique 

- Expand  parameter  database  for  ground  effects  modeling  for  low 
aspect  ratio  wings  (especially  LE  sweep  influence) 

• Computational  parametric  study  of  controlling 
factors 

- Unsteady  3-D,  to  study  LE  sweep  influence  on  DGE/SGE  differences 

- Steady-State  method,  check  static  ground  effect  levels 


During  this  test,  a list  of  items  that  need  to  be  reworked  was  generated.  Before 
this  DGE  cart  can  be  routinely  used  in  wind  tunnel  tests,  these  items  need  to  be 
fixed.  After  these  repairs  are  completed  (current  repair  plan  can  be  completed 
by  Feb/Mar  of  1999),  a follow-on  test  would  enable  us  to  check-out  the 
repairs.  After  a brief  check-out,  we  could  re-run  the  HSCT  planforms  already 
tested  to  confirm  the  current  findings  plus  run  the  seven  untested  models  to 
expand  the  ground  effects  database.  This  investigation  would  provide  the 
opportunity  to  gain  a better  understanding  of  the  differences  between  DGE  and 
SGE  test  techniques  as  well  as  the  knowledge  necessary  to  decide  when  each  is 
needed.  In  addition,  a computational  parametric  study  of  some  of  the  other 
potential  factors  would  help  the  HSR  community  gain  a better  understanding 
of  what  factors  need  to  be  controlled  in  ground  effects  testing. 
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Potential  Flow  Analysis 
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Task  33  4.3.2.4 


Potential  Flow  Analysis  of  Dynamic  Ground  Effect 
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Potential  Flow  Analysis  of  Dynamic  Ground  Effect 
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from  a camber  effect  because  of  the  flaps  being  down. 


F-15  Dynamic  Ground  Effect  Flight  Test  Results 

Figure  11.  Effect  of  approach  speed,  flap  setting,  and 
sink  rate  on  change  in  F-15  aerodynamic  coefficients 
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Lift  of  an  Impulsively  Started  Two-Dimensional  Flat  Plate  (Wagner  ’s  Function) 
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produced  by  virtual  mass  effects,  the  lower  leading  edge  suction  force  would  result  in  a reduced  amount  of  Polhamus  vortex 
lift  if  the  leading  edge  flow  is  separated. 
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C,c  (m)  = reference  chord  length 

H,h  (m)  = height  above  ground 

£2,(0  (rad/sec)  = frequency  of  oscillation 
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All  computations  were  performed  on  a 200MHz  PC.  The  code  was  compiled  using  the  Microsoft  Developer  Studio  FORTRAN  environment. 
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Validation  of  Time  Accurate  Vortex  Lattice  Method 
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Cropped  Delta  Wing  Lift  During  Constant  Angle-of- Attack  Descent  Towards  a Reflection  Plane 
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HSCT  High  Lift  Aerodynamics 

Cropped  Delta  Wing  Induced  Drag  During  Constant  Angle-of-Attack  Descent  Towards 

Reflection  Plane 
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HSCT  High  Lift  Aerodynamics 

Cropped  Delta  Wing  Induced  Drag  During  Constant  Angle-of- Attack  Descent  Towards  a 

Reflection  Plane 
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For  the  case  of  the  cropped  delta  wing,  all  lift  related  functions  proved  to  have  shapes  similar  to  a generalized  Wagner  function.  The  moment 
functions  were  nearly  independent  of  time  and  were  treated  as  constant  for  the  quick  preliminary  analyses  that  will  be  described  later. 


Unsteady  Pitching  Motion  of  Cropped  Delta  Wing 
Time  Accurate  Pitch  Stability  Derivatives 
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dCM/dq  time  dependent  pitch  damping 


HSCT  High  Lift  Aerodynamics 

Time  Dependent  Lift  and  Pitching  Moment  of  Cropped  Delta  Wing  Impulsively  Started  in  Free  Air 

at  10  Degrees  Angle-of-  Attack 
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HSCT  High  Lift  Aerodynamics 

Effect  of  Pitch  Rate  on  Dynamic  Lift  and  Pitching  Moment  of  Cropped  Delta  Wing  in 

Free  Air  at  Zero  Degrees  Angle-of-Attack 

Impulsive  Start.  Omegay=0.1,  Alpha  = 0 = constant. 
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The  figure  shows  that  the  airplane  will  have  traveled  nearly  4 chord  lengths  before  the  downward  lift  commanded  by  the  pilot  actually  begins 
develop.  For  a large  HSCT  traveling  at  150  knots  this  corresponds  to  a time  delay  of  about  1 second! 


G 

a> 

TD 

•O 

G £P 

C/5  G 

• 1*P 

£ ^ 

Gi 

O 0> 

~ Q 


O C/5 

pG  15 
u u 

s-»  <w 

Oh  g 


2353 


Quasi-Steady  Analysis  of  Pitch  Motion  Due  to  a 10  Degree  Sudden  Flap  Deflection  of 

a Neutrally  Stable  Cropped  Delta  Wing 
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TCA  Wing  Lift  Response  to  Sudden  Change  in  Ground  Height 
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8)  Tri-surface  configurations  exhibit  highly  complex  time  aerodynamic  responses  to  canard  deflections. 
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CPU  hours,  with  10  of  the  hours  required  to  obtain  the  initial  condition. 
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maKe  a aisconunuous  jump  to  tne  nigner  initial  lift  increment  values.  Instead  it  slowly  merged 
with  the  previous  analyses.  The  solution  continued  to  track  the  previous  solutions,  so  the 
analysis  was  stopped. 
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TCA  DGE  Results 
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Current  State  of  Affairs 
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Are  The  DGE  Characteristics  Tied  to  Viscous  Phenomena? 
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Recent  Results  in  the  Study  of 
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at  the  Biannual  Airframe  Technical  Review 
held  February  9-1 3, 1 998  at  the  Westin  Hotel, 
Los  Angeles  Airport 
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High  Lift  Technology  Development  (Task  33) 
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Programmed  Flaps  |_  Pressure  & Temp. 

Sensitive  Paints 
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GRIDTOOL;  the  initial  front  and  volume  grid  generation  in 
VGRID;  the  flow  solver  USM3D,  and  the  various  ways  used  to 
post-process  the  computational  results. 
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of  many  smaller  patches.  The  tetrahedra  are  background  sources 
that  can  be  used  to  control  the  size  and  distribution  of  the  tetrahedra 
during  the  volume  grid  generation. 
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body,  a result  of  the  background  sources  specified  in  the  initial  front 
development.  This  grid  distribution  served  as  a model  for  emulation 
when  other  wall  conditions  were  specified. 


I* 


a*5*?ss 


■fcftpaa 


&,v  ; : 

Bg&&iir~‘- 


Mfetoafi^i§g|gggkgj 


BfP 

BMi 


fe&’i  :^SSg9Bi 

EL  mi 

®jg&:  'ifl 
»g|^£35iS 


« s 

f-C  • *-H 

'4“*  X 

c 2 

- Oh 

> ^ 
> cd 


^ o_i 

c3  O 


CO 

C/3 

C+h 

s 

cd 

o 

o 

</3 

4— » 

0) 

3 

o 

G 

5-h 

o 

<d 

<D 

Sg 

too 

j> 

• i-H 

5-h 

<d 

> 

'o 

73 

c/3 

<D 

G 

£ 

5-H 

o 

O 

o 

O 

qG 

a 

<D 

5-h 

.2 

<u 

-G 

-G 

4— > 

73 

4— > 

1 

<D 

too 

g 

o 

1) 

G 

5-h 

* i-H 

G 

•i-H 

c 

cr 

CD 

G 

O 

5-h 

G 

C/3 

C/3 

i-H 

cd 

c 

<4-h 

O 

£ 

o 

•i-H 

4—* 

C/5 

♦ 

4— » 

cd 

3 

C/3 

<D 

H 

5-h 

<D 

4-J 

•i-H 

S-H 

o 

0) 

hC 

+-> 

o 

o 

H 

04  CO 

2482 


the  convergence  rate  to  this  number.  The  solution  was  accepted 
after  about  three  orders  of  magnitude  decrease  of  the  log  of  the 
residual. 
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coefficients  even  though  the  cell  sizes  did  vary.  Note  also  the  devas- 
tating result  of  removing  an  innocuous  side  source  on  the  number 
of  cells  in  the  volume. 


Initial  Grid  Study  Summary  17Jan98 
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of  the  wind  tunnel  box,  in  that  the  length  of  the  test  section  and  the 
position  of  the  model  were  modified  to  match  the  14-by-22  foot 
wind  tunnel  exactly. 
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aggravated  flow  in  that  area,  but  it  was  decided  to  go  with  the  origi 
nal  cell-counting  criterion  at  this  time.  A short  side  study  will  be 
done  to  look  at  the  effect  of  refining  the  grid. 


Real  Wind  Tdnnel”  Low-Position  Grid  Study  Summary  26Jan98 
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rotate  the  wind  tunnel  walls  down  8 degrees  to  present  a consistent 
flow  configuration.  Both  the  wind  tunnel  wall  and  the  plane  of  sym 
metry  surface  grids  are  shown  in  this  picture. 
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number  of  cells  remained  essentially  the  same  in  both  the  high  and 
low  cases  after  rotation,  so  these  grid  configurations  were  used  . 
Note  the  rotated  box  coordinates. 
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Wind  Tunnel  Data 
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point  is  slightly  in  error  for  the  USM3D  results,  and  it  is  anticipated 
that  movement  to  the  correct  location  will  eliminate  the  discrepan- 


TCA-1  WT  GE  Data,  LE=30,  TE=20 
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ment  were  examined  to  see  if  there  was  any  inconsistent  behavior 
of  the  data  at  low  model  heights  as  on  the  previous  page,  but  there 
did  not  seem  to  be  any. 
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Plans  for  the  Immediate  Future 
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HSCT  High  Lift  Aerodynamics 


2509 


ejr  - 2/10/98 


High  Lift  Technology  Development  (Task  33) 

% 

Increase  L/D,  Develop  Analysis/Design  Methodology 


KJ  jl. 

CD  O O < CL 


00 

<73 


E 

Cl 


2511 
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HSCT  High  Lift  Aerodynamics 
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influence.  It  is  important  that  the  traverser,  whose  purpose  is  to  measure  the 
flow  field  created  by  the  model,  not  alter  that  flow  field  and  give  inaccurate 
data.  Any  changes  to  the  flow  field  must  be  taken  into  account. 
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Possible  aerodynamic  influence 
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Operator  specified  traverser  path  and  boundary 
Portable 
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Single  point  attachment 


X W 0 ^ 

U-  o w = m *= 

O  ~ 0 £ ? O 

o — £ § ■*-  0 

■=0^20-0 

iElsf e 

EE  o £ 'S  •£ 

"O  < f?  O 05  g. 

!2  O ® ® :!  w 

3 ^ c £ 0 ^ 

M-  0 ^ C fc 

■=£30^0 

<2  +-  o -c  Q-  o 

0 g w c 5^  0 

S ® .s  S %>  S 

= S ® 5 s °- 

Q-  -C  "O  ZJ  Q_  0 

1 >*  5.  •-  E 0 

8 « ®’E»! 

0 ^ .£■  o ^ 

CO  Jr  L_  ° 0 0 

3 ^ "O  = ^ 

o 2 • 0 0 O 

.2  0 0 0 > O 

1 ! | o « ! 

.<2  ® o -g  « 

OT  £ E ® ® ■« 

^ ® t~5  ^ Cl 

« O £ f ® . 

10  § co  ^ 3 °-  -55 

£ .£  g c o £ o 

*o  "S  *2  Q-^-2  Q) 

0 ">  (A  0 8>‘“  C 

Q_  c =3  g £ Q-  0 

3 c 0 o m — 13 

a E 0 r e o) 

oSafill? 

-C  0 s-  O g > § 
h-  ID  O D £ $ O 


"O  C 
0 0 0 

0 £1 
-o  ^ «- 

C 0 o 

0-2.2  © 

5 o £ o c 

Zifa# 

® §l-,D§ 

32  U C 3 
£ 0 c/i  03  -Q 

E ^ c E w 

< c 0 0 

O o £ « £ 
p S?  ® -o 

S 2 E S'  S 
2 0 0 0 0 
£ Q-  o £ 0 

> - ? c°o 

® 8 <g  ■§  £ 

£ t:  00  0 O . 

^ 0 O O 0 0 

2  > 1 ® ® E 
® 8 >•  ® ■£-  8 

£ - 1 ° 3 s 
® to  § c | o 

5 g 0 C CD 
^ ^ -n  c c .2 

2 -2  0 2 "c  c 

£■05  s a ® 

nr  C > 2 c ° 

O 0 -rr  +3  ■—  Q. 

0 i_  m—  0 

•§  ® s ° £ ® 

m 9-  > 0 .0 

C 3 0 O 0 ° 
^ 0 O C _>,  0 

0 £ g 0 0 C 

3 X 2:  -c  c 0 

0 £ 0 X ^ E 

> o %p  21  V-  0 

1 -e  i b s 1 

Q-^EiS^o- 

^ ^ ‘o  P o P 


2522 


Rotating  the  lower  arm  would  have  required  re-paneling  the  base  of  the 
traverser.  To  avoid  this,  the  lower  arm  was  simply  kept  in  the  vertical  position. 
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Analysis  Configuration 
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where  this  problem  occurs,  mainly  over  the  leading  edge  flap  and  slightly  aft  of 
the  hingeline. 


2525 


of  leading  edge 
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horizontal 

Probe  tip  2.2  probe-diameters 
above  wing  surface  (0.6”) 
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uses  nigner-oraer  panel  method  to  solve  linearize 
potential  flow  boundary  value  problem 
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■Traverser  Effect  on  Pressures 
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Traverser  Effect  on  Forces  & Moments 
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High  Lift  Technology  Development  (Task  33) 
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Increasing  hingeline  radius  helps  the  flow  negotiate  the  upper  surface  curvature  at 
the  flap  break. 


HSCT  High  Lift  Aerodynamics 
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. Design  Study  Constraints 


0 

C 

o 

0 

1— 

0 

"O 

c 

o 

o 

0 

E 

-o 

c 

0 

0 

O 

0 


o 

-Q 

JD 

■D 

0 

C 

*0 
■*— » 
0 
C 

o 

o 

0 

0 


"D 

■4— » 

CO 

CD 


CO 

>,  E 

■D  0 

3 CO  C £ 

W | LL  w 


0 
-C  C 


_ 0 2 
o tr 


^ 0 c 

_ 0 O 0 

0 © 0)0 
■K  P .E  0 
W £ 0 E 
o E 0 


O -K 


0 


O)  E 
.E  o 


0 


E ° 
P co 


#o  8 S 

o1”  Jo 

£ g>  § S 

0'^-CQ- 
C 3 E -D 

O "D  $ 0 

.E  "P  CL  0 
^ 3 T3 


0 


■Q  0 


$ £ 0 c 
0 -D  g 8 

| 

*0  ■=  0 O 

> Z3  CL 

0 -Q  > "O 

s 5 2 ® 

•2  2 o c 

0 3 

> O 
0 O 0 


P 

3 

O 

£ 


0 

0 


0 
O)  0 


0 0 
■C  0 C 

S 0 .2> 

0 

$00 

0 E -o 

c o 


£ 

0 

> 


£ 

0 O 
C 


0 *-  P >,.E 
p o ^ c o 

I—  0 O)  0 Q. 


2560 


0 

E 

h- 


2561 


- potential  code  (A502) 


Model  Geometry 
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Inboard  Leading  Edges 
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launch  A502 


Sample  Gridding 
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1 . 35  degree  full  span  sealed  slats. 

2.  45  degree  full  span  sealed  slats. 

3.  35  degree  sealed  slats  outboard  and  30  degree  plain  flaps  inboard. 
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2570 


Neither  flap  configuration  is  over-deflected. 
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Neither  flap  configuration  is  over-deflected. 
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Both  configurations  have  a lower  surface  attachment  line;  there  was  some  concern  that 

commonly  deflected  full  span  flaps  / slats  would  result  in  an  over-deflection  of  the  inboard 
flaps. 
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about  55  counts. 
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The  photo  on  the  lower  right  shows  some  interference  between  the  new 
inboard  35  degree  sealed  slat  and  the  existing  outboard  35  degree  sealed  slat. 
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Navier-Stokes  Results  for 
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Lift  Coefficient  vs.  Drag  Coefficient 

Ref.  H,  Mach  0.24,  Alpha  10  Deg. 

Flaps  LE  Full-Span  30/30/30/30/30/30,  Partial-Span  0/0/30/30/30/30;  TE  10 
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Code  Calibration  Applied  to  the  TCA  High-lift 
Model  in  the  14x22  Wind  Tunnel 
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This  research  falls  under  the  Analysis/Design  Methodology  under  High  Lift  Development 
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The  outline  of  this  presentation  is  self-explanatory 
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HSR  High-Lift  TCA  Configuration 

Outline 
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The  objective  of  this  study  is  to  calibrate  a Navier-Stokes  code  for  the  TCA  (30/10)  baseline 


M 
i o 
CO)  C3 

e ts 

•j=i  c a 

2 'S 

5 s 

*£  'Eb 

c 

G 


G 

X 

G 

g 

CO 

W) 

0) 

X 

o 

cn 

4— > 

G 

X 

CD 

4— > 

o 

CD 

52 

CD 

X 

<D 

£ 

co 

& 

CD 

W> 

X 

D 

M 


G 

U 

CD 

> 

<D 

00 

U 

£ 

CO 

4— > 

3 

oo 

CD 

U 

g 

o 

* p*H 
4— > 

G 

4— » 

G 

cx 

s 

o 

o 

cd 

02 

H 


CO 

W) 

<D 


oo 

8 

G 

oo 

00 

CD 

U 

CX 

<D 

o 

X 

Uh 

u 

23 

00 

X 

C 

G 

rs 

CO 

4—* 

c 

<D 

S 


CD 

o 

U 

-O 


.5  73 

13  2 

<d 

f— H 4— » 

C * 

3 *0 


a, 

OO 

13 

t3 

a 

a, 

a 

o 


& 

o 

4> 

53 

4) 

"O 

2 

<D 


2 * 

2 & 
S «s 
So  53 
53  <u 
S bo 
o -o 

O 4) 


c3 

4— > 

g 

s 

•c 

CD 

CX 

X 

0) 


T3 

I 

CX 

S 

o 

0 

1 


X 

— CD 
52  co 

<5  3 


X) 

<d 

oo 

23 

8 

<D 

£ 


g 

x 


W) 

g 


G 

U 


D 


00 

CO 


00 
c 
O 
W)  co 


<d 
o 
a)  u 

X ,P 


c e 


s = 

G2  <5 


CO 

X 

<D 


X 

O 

o 

W) 

#s 

4> 


bQ  So 
• ~ 3 

-3  a> 
3 03 

g £ 


05 
05 

2 

G 

w)  EC 

c 


X 

g 

g 

bi) 

a 

o 

c 

D 

G 


X 

a 

G 

O 

Uf 

CX 

CX 

G 

X 

•c 


£ 

o 

g 


CD 
O 

G 

8 

CL> 

*— < 1+H 

w <d 

.g  06 

0 4) 
Cl,  jC 

1 4— > 

O 


G 

G 

O 


d 

X 


d 

4—* 

X3 

<D 


23 

co 

co 

G 

CO 

d 


CX 

6 

o 

o 

#s 

Vh 

<D 

> 

a) 

£ 

o 

X 

T3 

13 

T3 

O 

6 

d 


<u 

a 

ed 

0|— ( 

S 

CO 

Ui 

<D 

o 


<D 

CO 

u 

<D 

> 

2 

4) 


5 S 


4) 

JC 

4-J 

X3 

a> 

o 

c 

<u 

D 

55 

G 


cr  bo  c ^ 

d>  t *rn 


<D 

oo 

JC 

05 

4) 

s 

• #N 

Q 

CO 


2 Gi 


G 

6 


<u 

o 

G 

s? 

dj 

> 

G 

o 

o 

CD 


co  G 

•i-H  U 

JG  <D 


D 

a 

o 

G 


CX  §3 

^ g 

4> 

•c 

4) 
Oh 
X 
4) 


”0  -g 
4>  C 
C 


U 

05 

♦ 

>> 

XJ 

G 

w 

00 


X3 

CD 

CO 

G 

CD 

X) 

O 

o 

<D 

pG 

H 


o 

9 

G 

G 

u 

G 

bD 

52 

G 

O 

o 


<D 

O 

G 

(D 


co 

G 

CX 


G 

4—> 

X) 

o 

2 

4) 

£ 

co 

G 

O 

co 


O 

G 

co 

G 

£ 


a> 

o 

G 

<D 

X 


G 

G 

O 

52 

C 

W) 

* ^ 

CO 

4— > 

G 

G 

O 

s 

4— ► 

CO 

O 

CX 

<D 

X 

4— > 

X 

CD 


X 

O 

W) 

G 

•co 

G 


u 

W) 

W) 

X 

*4— ► 

CO 

• i-H 

X 

CD 

CD 

X 

4—* 

o 

4—) 

G 

X 

CD 

ts 

a> 

CO 

G 

CO 

G 


& 


G 

05  .2 


S & 


I 


CO 

<D 

o 

,G 

u 

CD 


X .(D 


<D 

4— * 

G 

G 

> 

X 

G 


CD 

C* 


X 
CD 
O 

« X 

■5  a> 

u 

O 53 
X CD 


o 

o 

8 

G 

co 

CO 

H 

CX 

X 

G 

G 


G 

<D 


<D 

f0> 

J+H 

Uh 

D 


G 

S 

CO 

D 

O 

G 

X 

O 

u 

CX 


X 

^G 

co 

2 

3 

CO 

CO 

8 

CX 

(D 

O 

.G 

U-i 

u 

G 

CO 

u 

<D 


S 3 
O o > 

g £ O 


0) 

g 

o 

4> 

bo 

4— » 

CO 

o 

CX 

D 

X 

o 

co 

CO 

2 

3 

CO 

CO 

2 

Oh 


05 

*o 

*c 

bo 


2696 


HSR  High-Lift  TCA  Configuration 

Objective  & Approach 
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tunnel  mount  effects. 
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HSR  High-Lift  TCA  Configuration 

Results 


0 

0 

© 


C/5 

© 


A -O 


0 

© 

P0 


© 

© 


bD 

<3 

73 

(N 


ia  *5 

o 

X) 

VO  00 
CM 

bO 

g 

• r-H 

G 

*s 

+-> 

g 
o 
o 


Vh 

bO 

M 

o 

0 

3 

1 


VO 


Vh 

a 

73 

<3 

a 

*3 

4-> 

X) 

o 

03 

g 

o 


r*  ^ a 


o 

00 


73 

<3 

03 

g 


t s. 


*r> 

VO 

<N 


73 

<3 

Vh 

• 

G 

CT 

<3 

Vh 


<3 


G 

O 

»i-H 

H-H 

G 

Vh 

(L) 


<3 

O 

”o 

<3 

03 


<N 


Vh 

O 

03 

Vh 

G 

O 

-G 

o 

<T5 

<3 

S 


a 

s 


<3 

73 

G 

4-> 

• 

G 

bX) 

G 

a 

3 

G 

73 

• t-H 

</3 

<3 

Vh 


<D 

03 

G 

<3 

Vh 

O 

<3 

73 

Vh 

<3 

73 

Vh 

O 


<3 

W>  _s 
G ^ 

Vh  , 

<3  1 

> 

<1  <N 


I I I I 


2699 


3 

X 

H 


3 

7D 

S-N 

3 

a 

X 

W) 

o 

X 

cd 

C/3 

CL 

IS 

m 

■w 

73 

c 

a 

* ^ 

cd 

c 

£ 

C/3 

3 

o 

o 

X 

,cd 

C/3 

C/3 

(3 

• ^ 

-4-J 

3 

c 

c 

cd 

OX) 

"CL 

c 

Ui 

4— > 

3 

X 

o 

4— » 

cd 

£ 

£ 

s 

4— » 

c 

• ^4 

C/1 

o 

3 

Cl 

X 

*-* 

6 

4—* 

c 

1 

4-4 

* *-< 

c 

73 

3 

4— > 

* *— < 
o 

CL 

O 

m 

U 

U7 

73 

c/T 

’C 

CJO 

X 

CJ 

0) 

o 

3 

X 

i-H 

x 

CN 

p 

L— i 

c ^ 

o 

o 

73 

3 

C/3 

o 

O 

co 

CL 

S 

< 

o 

CJ 

3 

H 

C/3 

4> 

3 

sz 

*c 

H 

M 

so 


> 

o 


X 

6 

d 

i 

u 

03 

73 

0) 

cn 

P 


c/5 

3 


73 

c 

P 

O 

X) 


c/5 

•4— > 

c 

'o 

CL 


2700 


TCA  (30/10)  Baseline  Grid  With  No  Post 
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Spanwise  Pressure  Distributions  for  TCA  (30/1 0) 
With  No  Post-Mount 

M = 0.243,  a = 8°,  Re  - 0 mil 
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The  post  effects  were  considered  next;  Allen  Chen  from  Boeing  (using  the  panel  code 
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HSR  High-Lift  TCA  Configuration 

Results 
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A global  view  of  the  TCA  (30/10)  configuration  mounted  on  the  post  is  shown.  The  post 
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A close-up  view  of  the  model  support,  which  interfaces  with  the  lower  wing  surface.  The 
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Method  of  Attaching  Post-Mount  to  TCA  Model 
Using  Overset  Grids 
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A more  global  view  of  the  TCA  (30/10)  surface  grid  and  post  is  shown  on  the  left  of 
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The  pressure  contours  are  plotted  in  the  symmetry  plane  and  on  the  lower  configuration 
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Pressure  Contours  Plotted  in  the  Overset  Grid  Regions 
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TCA  Baseline  (30/10)  Configuration 

Lower  Surface  C Contours  With  and  Without  Post-Mount 
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The  chordwise  pressure  distributions  with  and  without  the  post-mount  are  compared  at 
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Chordwise  Pressure  Distributions  for  TCA  (30/10) 
With  and  Without  Post-Mount 

M = 0.243,  a = 8°,  Re  = 8 mil 
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Spanwise  Pressure  Distributions  for  TCA  (30/1 0) 
With  and  Without  Post-Mount 

M = 0.243,  a = 8°,  Re  = 8 mil 
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HSCT  High  Lift  Aerodynamics 
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CL  Comparisons  With  and  Without  the  Post-Mount 

TCA  (30/1 0)  Baseline,  M = 0.243,  Re  = 8 mil. 
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Comparisons  With  and  Without  the  Post-Mount 

TCA  (30/10)  Baseline,  M = 0.243,  Re  = 8 mil. 
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Drag  Polar  Comparisons  With  and  Without  the  Post-Mount 

TCA  (30/10)  Baseline,  M = 0.243,  Re  = 8 mil. 
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Comparisons  With  and  Without  the  Post-Mount 

TCA  (30/10)  Baseline,  M = 0.243,  Re  = 8 mil. 
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HSR  High-Lift  TCA 

Summary  & Conclusions 
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Aerodynamic  Analysis  of  TCA 
Wing/Body/Nacelle  High  Lift  Configurations 


Xuetong  Fan 
Paul  Hickey 
ASE  Technologies,  Inc. 

High  Speed  Research  Program 
Airframe  Technical  Review 
Westin  Hotel,  Los  Angeles  Airport 
February  9-13,  1998 


This  presentation  includes  the  work  completed  at  ASE  Technologies,  Inc. 
in  the  CFD  analysis  of  HSCT  TCA  Wing/Body/Nacelle  High  Lift  configurations. 
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High  Lift  Technology  Development  (Task  33) 

Increase  L/D,  Develop  Analysis/Design  Methodology 


In  the  High  Lift  Technology  Development  program  (Task  33), 
this  work  falls  in  the  category  of  Navier-Stokes  Methods 
under  Analytical  Methods.  We  hereby  acknowledge  the 
support  and  help  from  Roger  Clark  and  David  Yeh  of  Boeing 
Long  Beach. 
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Objectives 


• Develop  CFD  Models  for  aerodynamic  analysis  of  HSCT  TCA 

Wing/Body/Nacelle  (WBN)  high  lift  configurations 

- Construct  multi-block  CFD  grids  to  include  nacelles  and  diverters 
with  and  without  deflected  LE/TE  flaps 

- Obtain  converged  CFL3D  solutions  for  the  TCA  WBN  models 

• Evaluate  the  effect  of  nacelle  installation  on  the  aerodynamic 

performance  of  HSCT  TCA  high  lift  configurations 

- Identify  and  analyze  important  flow  characteristics  due  to  nacelle 
installation  to  support  Propulsion  Airframe  Integration  (PAI)  effort 

- Provide  flow  and  performance  data  related  to  nacelle  effect  to 
supplement  wind  tunnel  test 


The  objectives  of  this  work  are  two  fold.  The  first  objective  is  to  develop 
efficient  CFD  modeling  procedures  for  the  TCA  high  lift  configurations  with 
nacelle  installation.  The  second  objective  is  to  evaluate  the  effect  of  nacelle 
installation  on  the  aerodynamic  performance  of  the  TCA  high  lift  configurations. 

To  achieve  the  first  objective,  we  will  build  multi-block  CFD  grids  to  include 
nacelles  and  diverters  in  the  TCA  high  lift  configurations  with  and  without 
deflected  LE/TE  flaps.  And  then  we  will  use  CFL3D  to  obtain  fully  converged 
turbulent  solutions  for  the  TCA  W/B/N  models. 

For  the  second  objective,  we  will,  from  the  CFD  solutions,  identify  and  analyze 
important  flow  characteristics  due  to  nacelle  installation  to  support  Propulsion 
Airframe  Integration.  We  will  also  provide  flow  and  performance  data  for  the 
TCA  W/B/N  configurations  to  supplement  wind  tunnel  test. 
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Outline 

• TCA  0/0  WBN  Model 

- Model  description 

- Convergence  history 

- Comparison  with  wing/body  solution  and  test  data 

- Flow  characteristics 

• TCA  30/10  WBN  Model 

- Model  description 

- Convergence  history 

- Comparison  with  wing/body  solution  and  test  data 

- Flap  loading  analysis 

- Flow  characteristics 

• Summary 


We  completed  two  CFD  models  for  the  TCA  W/B/N  configurations:  TCA  0/0 
and  TCA  30/10.  In  this  presentation,  we  will  describe  the  multi-block  CFD 
models  for  the  two  W/B/N  configurations  and  show  the  CFL3D  solutions. 

We  will  compare  our  CFL3D  solutions  for  W/B/N  with  the  respective  CFL3D 
solutions  for  W/B  obtained  at  Boeing  Long  Beach  for  the  same  flap  settings. 
We  will  also  compare  the  CFD  solutions  with  the  available  test  data  to  illustrate 
nacelle  effect.  We  will  use  flow  visualization  to  show  the  important  flow 
characteristics,  especially  in  the  vicinity  of  nacelles  and  diverters.  For  the 
TCA  30/10  W/B/N  configuration,  we  will  examine  the  span  loading  of  the 
deflected  TE  flaps.  A brief  summary  of  this  work  is  included  at  the  end  of 
this  presentation. 
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TCA  0/0  WBN  Model  Description 


Based  on  single-block 
TCA  0/0  W/B  grid 
provided  by  David  Yeh 

Nacelle/Diverter  (N/D)  grid 
face-matched  with  W/B  grid 

Mostly  1 — to— 1 point-matched 
between  N/D  blocks 

Model  size:  20  blocks 
4 million  grid  points 

Memory  requirement:  160MW 

Minimum  viscous 
spacing:  0.002  in 


The  construction  of  TCA  0/0  W/B/N  grid  is  based  on  the  TCA  0/0  W/B  grid 
provided  by  David  Yeh.  A portion  of  the  W/B  grid  is  removed  to  make  room 
for  nacelles  and  diverters  (N/D).  The  N/D  blocks  interface  with  the  W/B  grid 
in  a face-matched  manner  using  the  RONNIE  pre-processor  in  CFL3D. 
Between  the  N/D  blocks,  1 — to— 1 point-match  is  used  as  much  as  possible. 
The  final  model  for  TCA  0/0  W/B/N  has  20  blocks  with  4 million  grid  points. 
It  requires  1 60  MW  of  memory  on  Cray  C-90  to  run  the  CFL3D  flow  solver 
for  this  model.  The  minimum  viscous  spacing  in  the  W/B/N  model  is  0.002 
inches  to  be  consistent  with  the  W/B  grid. 
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TCA  0/0  WBN  Model 


With  the  minimum  viscous  spacing  of  0.002  inches,  the  y+  values  for  the  first 
grid  points  off  model  solid  surfaces  are  in  the  order  of  1 .0,  which  is  appropriate 
for  the  Baldwin-Lomax  turbulence  model.  The  flow  condition  for  the  y+  results 
are:  Mach  No.  = 0.3,  Re  = 8 million,  and  Angle  of  Attack  = 10  degrees. 
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CFL3D  Convergence  History 

(TCA  0/0,  Mach  #=0.3,  Re=8x106,  AOA=1 0 degrees) 


Coarse  Level  Fine  Level 


• 2-level  grid  sequencing;  2-level  multi-grid  on  fine  level 

• Total  CPU  Usage:  42  hours  on  Cray  C-90 

(9  hours  on  coarse  level  and  33  hours  on  fine  level) 


The  convergence  history  of  the  CFL3D  solution  for  the  TCA  0/0  W/B/N  model 
is  shown  here.  Two-level  grid  sequencing  was  used  and  two-level  multi-grid 
was  applied  on  the  fine  grid  level.  Extensive  iterations  were  performed  on  the 
coarse  grid  level  for  debugging  purposes,  which  used  9 hours  of  CPU  time  on 
Cray  C-90.  After  500  iterations  on  the  fine  level,  we  adopted  David  Yeh’s 
modification  to  the  Baldwin-Lomax  turbulence  model  with  Degani-Schiff  option. 
It  has  been  proven  that  Yeh’s  modification  can  better  simulate  the  vortical  flow 
above  the  wing  upper  surface.  For  simplicity,  we  only  used  his  modification  in 
the  block  around  the  leading  edge  and  above  the  wing  upper  surface  in  the 
W/B/N  model.  This  modification  caused  the  predicted  lift  coefficient  to  increase 
gradually,  which  is  why  extensive  (and  expensive)  iterations  were  run  on  the 
fine  grid  level.  The  1800  iterations  on  the  fine  grid  level  used  33  CPU  hours 
on  Cray  C-90. 
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CFL3D  Convergence  History 

(TCA  0/0,  Mach  #=0.3,  Re=8x106,  AOA=10  degrees) 


The  convergence  history  for  CL  and  CD  for  the  TCA  0/0  W/B/N  model 
is  shown  here.  For  the  last  200  iterations,  the  peak-to-valley  variations 
in  CL  and  CD  are  within  2%  of  their  mean  values,  which  are  comparable 
to  the  level  of  convergence  achieved  in  the  W/B  model.  Since  our  W/B/N 
model  is  based  on  the  W/B  grid,  it  is  unlikely  for  the  W/B/N  solution  to 
converge  much  better  than  the  W/B  solution. 
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Comparison  with  W/B  Solution  and  Wind  Tunnel  Data 

(TCA  0/0,  Mach  #=0.3,  Re=8x106,  AOA=10  degrees) 

Lift  Coefficient  (Cl) 


CL 

W/B 

W/B/N 

ACL 

CFL3D 

0.3850 

0.4116 

0.0266 

Experiment 

0.3749 

0.371 1 

-0.0038 

Drag  Coefficient  (Cp) 


Cd 

W/B 

W/B/N 

ACd 

CFL3D 

0.0599 

0.0634 

0.0035 

Experiment 

0.0606 

0.0616 

0.0010 

The  integrated  coefficients,  CL  and  CD,  for  the  TCA  0/0  W/B/N  model  are 
shown  in  the  above  tables,  together  with  those  from  the  W/B  solution  and 
wind  tunnel  experiments.  Considerable  discrepancies  are  found  between 
the  CFD  solutions  and  the  experimental  data.  Further  and  more  detailed 
investigation  is  necessary  to  understand  and  possibly  eliminate  the 
discrepancies.  It  is  noted  that  the  Cd  values  from  both  the  W/B/N  CFD 
solution  and  the  test  data  have  been  corrected  for  the  nacelle  interior 
friction  forces  and  the  nacelle  base  pressure  forces. 
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Comparison  with  W/B  Solution  and  Wind  Tunnel  Data 

(TCA  0/0,  Mach  #=0.3,  Re=8x106,  AOA=10  degrees) 


This  figure  shows  the  relative  locations  of  the  lift  coefficients  on  the  lift  curve. 

It  is  noted  that  for  the  angles  of  attack  (AoA)  up  to  16  degrees,  wind  tunnel  test 
shows  decreases  in  C^,  whereas  the  CFD  solution  predicts  an  opposite  trend 
at  AoA  of  1 0 degrees. 
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This  slide  shows  the  pressure  coefficient  contour  on  the  wing  upper  surfaces 
for  the  W/B  and  W/B/N  CFD  models.  With  nacelles  installed,  the  CFD  solution 
shows  similar  leading  edge  vortices  but  with  enhanced  strength,  compared  to 
the  W/B  solution.  This  enhancement  resulted  in  the  increment  in  the  predicted 
lift  coefficient  for  the  TCA  0/0  WBN  configuration. 
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On  the  wing  lower  surface,  the  pressure  coefficient  distributions  and  the 
limiting  streamlines  are  very  similar  between  W/B  and  W/B/N  solutions, 
except  for  around  the  diverters.  For  the  W/B/N  case,  the  limiting  streamlines 
indicate  possible  flow  separation  near  the  leading  edges  of  the  diverters. 
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TCA  0/0  WBN  CFL3D  Solutions: 

(Mach  #=0.3,  Re=8x106,  AOA=10  degrees) 

Normalized  total  pressure  downstream  of  nacelle  inlets 


Inboard 


Pt 


— 1.05 


Outboard 


This  figure  shows  the  total  pressure  contours  in  the  nacelle  interior  cross 
sections  just  downstream  of  the  nacelle  leading  edge.  Local  flow  separation 
is  apparent  at  the  entrance  to  the  outboard  nacelle  which  may  impact 
the  outboard  engine  performance. 
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TCA  0/0  WBN  CFL3D  Solutions: 

(Mach  #=0.3,  Re=8x106,  AOA=10  degrees) 
Limiting  streamlines  around  nacelle  leading  edge 


Examination  of  the  limiting  streamlines  near  the  leading  edge  of  the 
outboard  nacelle  further  confirms  the  local  flow  separation. 
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The  W/B/N  model  for  the  TCA  30/1 0 configuration  is  based  on  the  TCA  30/1 0 
W/B  grid  provided  by  David  Yeh.  To  improve  the  interfacing  quality  between 
the  W/B  grid  and  the  N/D  grid,  we  used  buffer  zones  between  them  at  spanwise 
locations.  These  buffer  zones  1 — to— 1 point-match  with  the  W/B  grid  and 
face-match  with  the  N/D  grid  on  flat  surfaces  for  higher  interpolation  accuracy. 
Between  the  N/D  blocks,  mostly  1 — to— 1 point-match  is  used. 

The  major  technical  challenge  for  the  TCA  30/10  W/B/N  configuration  is  the 
modeling  of  the  deflected  TE  flaps.  With  nacelles  installed,  the  "web"  approach 
used  by  David  Yeh  for  the  W/B  model  is  not  feasible.  In  addition,  due  to  the 
close  proximity  of  the  deflected  TE  flaps  to  the  nacelles  and  diverters,  it  is 
difficult  to  obtain  high  quality  grid  that  can  accurately  simulate  the  potentially 
complicated  flows  in  that  region. 
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TCA  30/10  WBN  Model 


Our  approach  to  the  problem  is  to  use  three  wedges  for  each  region  above 
the  upper  surfaces  of  the  deflected  flaps,  as  shown  in  the  above  figure. 
These  wedges  fill  in  the  space  created  by  TE  flap  deflection  and  at  the  top 
flush  with  the  wing  upper  surface.  The  side  surfaces  of  the  undeflected  wing 
segments  and  the  openings  between  the  wing  and  the  deflected  TE  flaps 
are  accurately  modeled  with  this  approach. 
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TCA  30/10  WBN  Model 

Modeling  the  deflected  TE  flaps  with  wedges 


Underneath  the  wing,  another  three  wedges  are  used  for  each  region  between 
the  deflected  TE  flaps.  The  bottom  surface  of  these  wedges  are  flush  with 
the  lower  surfaces  of  the  deflected  TE  flaps.  Through  the  openings  created  by 
the  TE  flap  deflections,  1 — to — 1 point-match  is  established  between  the  upper 
and  lower  wedges. 
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TCA  30/10  WBN  Model 

• Model  size:  50  blocks  (including  21  wedges);  6 million  grid  points 

• Memory  requirement:  255MW 


This  figure  shows  the  actual  wedge  grid  in  one  of  the  deflected  TE  flap  regions. 
It  is  noted  that  one  side  for  each  set  of  the  lower  wedges  will  lie  on  the  side 
surfaces  of  the  nacelles  and  diverters.  The  final  grid  for  the  TCA  30/10  W/B/N 
configuration  consists  of  50  blocks  including  21  wedges  and  a total  of  about 
6 million  grid  points.  It  requires  255  MW  of  memory  on  Cray  C-90  to  run  the 
CFL3D  solver  for  this  model. 


2753 


A S E Technologies , Inc , 


CFL3D  Convergence  History 

(TCA  30/10,  Mach  #=0.3,  Re=8x106,  AOA=10  degrees) 


iterations 

1^ >1^-4 >1 

Coarse  Level  Intermediate  Fine  Level 
Level 

• 3-level  grid  sequencing 

• 2-level  multi-grid  on  intermediate  and  fine  levels 

• CPU  Usage:  54  hours  on  Cray  C-90  (2  hours  on  coarse  level; 

2 hours  on  intermediate  level;  50  hours  on  fine  level) 


With  the  improved  accuracy  for  the  patched  grid  interfaces  in  the  TCA  30/1 0 
W/B/N  model,  we  are  able  to  run  3-level  grid  sequencing  with  2-level  multi-grid 
on  the  intermediate  and  fine  grid  levels.  We  completed  3000  iterations  on  the 
coarse  grid  level  and  600  iterations  on  the  intermediate  level  which  used 
2 CPU  hours  each  on  Cray  C-90.  After  switching  to  fine  grid  level  after  3600 
iterations,  however,  fluctuations  started  to  appear  in  the  residual  history. 

Various  means  were  tried  to  elliminate  the  fluctuations  with  no  success. 
Examination  of  the  flow  field  revealed  that  local  flow  fluctuations  exist  only  in 
the  wedge  regions.  After  6300  total  iterations,  we  adopted  David  Yeh’s 
modification  to  the  turbulence  model  and  completed  another  300  iterations. 

All  the  fine  level  iterations  used  50  CPU  hours  on  Cray  C-90. 
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CFL3D  Convergence  History 

(TCA  30/10,  Mach  #=0.3,  Re=8x106,  AOA=10  degrees) 


Iterations 


Iterations 


This  figure  shows  the  convergence  history  for  the  integrated  lift  and  drag 
coefficients  for  TCA  30/1 0 W/B/N  model.  The  local  fluctuations  in  the  flow 
field  near  the  deflected  TE  flaps  have  no  apparent  effect  on  the  integrated 
coefficients.  In  addition,  Yeh’s  modification  to  the  turbulence  model  is  not 
affecting  the  solution  due  to  the  weakness  of  the  LE  vortices  for  the  TCA 
30/1 0 configurations.  For  the  last  200  iterations,  the  peak-to-valley  variations 
in  CL  and  CD  are  within  2%  of  their  mean  values  respectively,  which  are 
comparable  to  the  level  of  convergence  achieved  in  the  TCA  30/10  W/B 
CFD  model. 
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Comparison  with  W/B  Solution  and  Wind  Tunnel  Data 

(TCA  30/10,  Mach  #=0.3,  Re=8x106,  AOA=10  degrees) 

Lift  Coefficient  (Cl) 


CL 

W/B 

W/B/N 

ACL 

CFL3D 

0.4421 

0.4635 

0.0214 

Experiment 

0.4697 

0.0437 

Drag  Coefficient  (Cp) 


Cd 

W/B 

W/B/N 

ACq 

CFL3D 

0.0501 

0.0571 

0.0070 

Experiment 

0.0505 

0.0602 

0.0097 

The  integrated  coefficients,  CL  and  CD,  for  the  TCA  30/10  W/B/N  model  are 
shown  in  the  above  tables,  together  with  those  from  the  W/B  solution  and 
wind  tunnel  experiments.  Discrepancies  are  found  between  the  CFD  solutions 
and  the  experimental  data  in  terms  of  lift  and  drag  increments  due  to  nacelle 
installation.  Further  and  more  detailed  investigation  is  necessary  to  understand 
and  possibly  eliminate  these  discrepancies.  It  is  noted  that  the  Cq  values  from 
both  the  W/B/N  CFD  solution  and  the  test  data  have  been  corrected  for  the 
nacelle  interior  friction  forces  and  nacelle  base  pressure  forces. 


2756 


A SE  Technologies , Inc . 


Comparison  with  W/B  Solution  and  Wind  Tunnel  Data 

(TCA  30/10,  Mach  #=0.3,  Re=8x106,  AOA=1 0 degrees) 


a 


This  figure  shows  the  relative  locations  of  the  lift  coefficients  on  the  lift  curve. 
For  TCA  30/10  configurations,  both  CFD  analysis  and  wind  tunnel  experiment 
show  an  increase  in  CL  due  to  nacelle  installations,  although  the  CFD  solution 
predicts  a much  smaller  increment  in  CL  than  test  data  at  AoA  of  1 0 degrees. 
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Comparison  with  W/B  Solution  and  Wind  Tunnel  Data 

(TCA  30/10,  Mach  #=0.3,  Re=8x106,  AOA=10  degrees) 


This  figure  shows  the  relative  locations  of  the  lift  and  drag  coefficients  on 
the  drag  polar  for  the  TCA  30/1 0 W/B  and  W/B/N  configurations. 
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Comparison  with  W/B  Solution 

(TCA  30/10,  Mach  #=0.3,  Re=8x106,  AOA=10  degrees^ 
Cp  on  wing  upper  surface 

Wing/Body 


Wing/Body/Nacelle 


The  pressure  coefficient  contour  on  the  wing  upper  surface  for  the  TCA 
30/10  W/B/N  model  is  similar  to  that  for  the  W/B  model  due  to  the  weakness 
in  the  vortices  around  the  deflected  leading  edge.  The  nacelles  have  very 
little  effect  on  the  overall  pressure  distribution  on  the  wing  upper  surface. 
Local  "hot  spots"  can  be  seen  near  the  deflected  TE  flaps  which  are 
responsible  for  the  residual  fluctuations  in  the  convergence  history. 
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On  the  wing  lower  surface,  the  pressure  coefficient  distributions  and  the 
limiting  streamlines  are  very  similar  between  the  W/B  and  W/B/N  solutions, 
except  for  around  the  diverters.  For  the  W/B/N  case,  the  limiting  streamlines 
indicate  possible  flow  separation  near  the  leading  edge  of  the  diverters.  Local 
flow  acceleration  can  be  seen  between  the  inboard  and  the  outboard  nacelles. 
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TCA  30/10  WBN  CFL3D  Solutions: 

(Mach  #=0.3,  Re=8x1 06,  AOA=10  degrees) 
Normalized  total  pressure  downstream  of  nacelle  inlets 


Inboard 


pt 


Outboard 


This  figure  shows  the  total  pressure  contours  in  the  nacelle  interior  cross 
sections  downstream  of  the  nacelle  leading  edge.  Local  flow  separation 
is  apparent  at  the  entrance  to  the  outboard  nacelle  which  may  impact 
the  outboard  engine  performance.  Flow  separation  also  exists,  though 
not  severe,  at  the  entrance  to  the  inboard  nacelle. 


2761 


c C,  /c 


ASE  Technologies , Inc . 


Comparison  with  W/B  Solution:  Flap  Span  Loading 

(TCA  30/10,  Mach  #=0.3,  Re=8x106,  AOA=10  degrees) 


Wing/Body Wing/Body/Nacelle 


This  slide  shows  the  effect  of  nacelles  on  the  span  loading  of  the  deflected 
TE  flaps.  It  is  expected  that  the  nacelles  will  work  like  end-plates  to  the 
wing  or  flap  segments,  which  will  usually  cause  an  increase  in  span  loading. 
This  end-plating  effect  is  apparently  seen  for  the  inner  and  outer  flaps 
but  not  for  the  middle  flap. 
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TCA  30/10  WBN  CFL3D  Solutions: 

(Mach  #=0.3,  Re=8x106,  AOA=10  degrees) 
Cp  on  lower  surfaces  of  deflected  middle  flaps 


Wing/Body 


Wing/Body/Nacelle 


Further  examination  of  the  pressure  distribution  on  the  lower  surface 
of  the  deflected  middle  flap  shows  relatively  lower  overall  pressure 
for  W/B/N  case  compared  to  the  W/B  case.  This  is  the  result  of 
local  flow  acceleration  between  the  inboard  and  outboard  nacelles. 
The  flow  acceleration  and  end-plating  effect  tend  to  offset  each  other 
in  the  span  loading  on  the  middle  flap. 
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Summary 

• Developed  CFD  Models  for  TCA  Wing/Body/Nacelle  (W/B/N) 
high  lift  configurations  including  0/0  and  flaps  30/1 0. 

• Obtained  converged  CFL3D  solutions  for  the  TCA  W/B/N  models 
at  Mach  No.=0.3,  Re=8x106,  AOA=10  degrees. 

• CFD  solutions  provide  insight  into  the  effect  of  nacelle  installation 
on  the  aerodynamic  performance  of  TCA  high  lift  configurations. 

- Impact  on  overall  flow  field  and  aerodynamic  parameters 

- Local  flow  separation  on  nacelles/diverters 

- Span  loading  and  flow  characteristics  near  deflected  TE  flaps 

• Comparison  between  CFD  results  and  wind  tunnel  data  reveals 
discrepancies  in  aerodynamic  coefficients.  Further  investigation 
is  warranted. 


In  summary,  we  developed  two  CFD  models  and  obtained  converged 
CFL3D  solutions  for  the  TCA  Wing/Body/Nacelle  high  lift  configurations 
including  flaps  0/0  and  30/10.  CFD  solutions  make  it  possible  to  examine 
in  detail  and  understand  in  depth  the  effect  of  nacelle  installation  on  the 
aerodynamic  performance  of  HSCT  TCA  high  lift  configurations. 
Preliminary  comparison  with  wind  tunnel  test  data  reveals  discrepancies 
in  the  integrated  aerodynamic  coefficients.  We  will  perform  further 
investigation  into  the  CFD  solution  as  well  as  the  test  data. 
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High  Lift  Technology  Development  (Task  33) 

Increase  UD,  Develop  Analysis/Design  Methodology 
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Outline 
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Approach 
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W/B/C  Modeling  Using  Overset  Approach 
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W/B/C  Modeling  Using  Overset  Approach 
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Automated  Canard  Modeling  Procedure 
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Automated  Canard  Modeling  Procedure 
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Aft  elevator  modeling 
Forebody/wiping  surface 


Automated  Procedure  for  Canard  Surface  Integration 
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Automated  Procedure  for  Canard  Surface  Integration 

(T)  Input  - Surface  Grids  (f)  Surface  Grid 
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Create  canard  grid 
(on  exposed  area) 


Surface  Geometry  Definition  for  Modeling  Canard  Surface 
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Surface  Geometry  Definition  for  Modeling  Canard  Surfaces 

HSCT  Aerodynamics,  Long  Beach 
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Canard  Manipulation  for  Sensitivity  Studies 
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Canard  Manipulation  for  Sensitivity  Studies 

HSCT  Aerodynamics,  Long  Beach 
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Canard  Translation  Process 
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Canard  Dihedral/Anhedral  Setting  Modification 
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Canard  Dihedral/Anhedral  Setting  Modification 

■ ■■  i.  n i umm  HSCT  Aerodynamics,  Long  Beach 
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Canard  Hingewise  Deflection 
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Canard  Hingewise  Deflection 
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Numerical  Procedure  for  Trim  Curve  Calculation 
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Numerical  Procedure  for  Trim  Curve  Calculation 

HSCT  Aerodynamics,  Long  Beach 
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Surface  Grid  Construction  Procedure 
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Surface  Grid  Construction  Procedure 
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Sectional  Volume  Grid  Generation 
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Sectional  Volume  Grid  Generation 
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Forebody/Canard  Volume  Grid 
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Numerical  Integration  Procedure  for  W/B/C  Configurations 
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Create  Overset  Connectivity  Using  MAGGIE 
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TCA  Wing/Body/Canard  Geometry 

TCA  Wing/Body  with  PTC  (v2)  Canard,  i = 0.0° 
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Flow  Solution  for  TCA  W/B/C  Configuration 
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Flow  Solution  for  TCA  W/B/C  Configuration 

HSCT  Aerodynamics,  Long  Beach 
M=0.3,  AoA=10  deg,  Re=8  million 
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Flow  Solution  for  TCA  W/B/C  (30/10)  Configuration 
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Canard  Influence  on  the  Wing  Upper  Surface  Pressure 
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Canard  Influence  on  the  Wing  Limiting  Streamlines 
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Effect  of  Canard  on  Wing  Surface  Pressure 
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Effect  of  Canard  on  TCA  (30/10)  Wing  Separation 
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Effect  of  Canard  on  TCA  (30/10)  Wing  Separation 

HSCT  Aerodynamics,  Long  Beach 

M=0.3,  AoA=10  deg,  Re=8  million 


Canard  configuration  for  TCA-4  High  Lift  Test 
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HSCT  Aerodynamics,  Long  Beach 
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• Influence  of  PTC  canard  on  the  wing  should  be  less  than  ACC  canard 
due  to  its  size  and  position 


High  Lift  CFD  Activities  for  Canard  Integration 
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- Canard  planform 

- Canard  position 
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Comparison  of  CFL3D  Solutions  Using 
Alternative  Grid  Interfacing  Schemes 


Xuetong  Fan 
Paul  Hickey 
ASE  Technologies,  Inc. 

High  Speed  Research  Program 
Airframe  Technical  Review 
Westin  Hotel,  Los  Angeles  Airport 
February  9-13, 1998 


This  presentation  documents  the  recent  work  at  ASE  Technologies,  Inc. 
in  applying  the  overset  grid  method  to  the  CFD  modeling  of  HSCT  TCA 
Wing/Body/Nacelle  High  Lift  configurations.  The  MAGGIE  preprocessor 
in  CFL3D  is  used  in  this  work  for  the  grid  overlapping  interpolations. 


2831 


High  Lift  Technology  Development  (Task  33) 

Increase  UD,  Develop  Analysis/Design  Methodology 


In  the  High  Lift  Technology  Development  program  (Task  33), 
this  work  falls  in  the  category  of  Navier-Stokes  Methods  under 
Analytical  Methods.  We  hereby  acknowledge  the  support  and 
help  from  Roger  Clark  and  David  Yeh  of  Boeing  Long  Beach. 
We’d  also  like  to  thank  Chung-Jin  Woan  of  Boeing  North 
American,  Inc.  for  the  fruitful  discussions  on  the  overset  grid 
method. 
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Objectives 

• Compare  CFL3D  solutions  for  TCA  Wing/Body/Nacelle  (WBN) 
high  lift  configurations  using  different  grid  interfacing  schemes 

- Face-matched  grid  (RONNIE) 

- Overlapped  grid  (MAGGIE) 

• Evaluate  grid  interfacing  schemes  for  more  efficient  CFD  modeling 
of  TCA  WBN  configurations  with  and  without  deflected  LE/TE  flaps 

- Grid  generation  effort 

- Grid  interface  quality 

- Computer  resources 

- Technical  issues 


The  objectives  of  this  work  are  twofold.  The  first  objective  is  to  compare 
and  cross-examine  the  CFL3D  solutions  for  the  TCA  W/B/N  models  using 
two  different  grid  interfacing  methods:  face-matching  (or  patching)  with 
RONNIE  and  overlapping  with  MAGGIE.  The  second  objective  is  to  evaluate 
these  two  methods  and  determine  which  one  is  better  suited  for  the  CFD 
modeling  of  the  TCA  W/B/N  configurations.  The  grid  interfacing  method  will 
be  evaluated  in  terms  of  grid  generation  effort,  block  interface  quality,  and 
computer  resources.  Special  attention  is  paid  to  the  potential  technical 
difficulties  involved,  especially  in  the  case  of  deflected  TE  flaps  for  TCA  high 
lift  configurations. 
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Outline 

• Background 

• Overlapped  grid  model  for  TCA  0/0  WBN  High  Lift  configuration 

• CFL3D  solution  for  TCA  0/0  WBN  model  with  overlapped  grid 

- Computer  resources 

- Comparison  with  CFL3D  solution  using  face-matched  grid 

• Comparison  of  the  two  schemes  in  TCA  WBN  applications 

• Technical  issues  with  overlapped  grid 

- Interpolation  near  solid  surfaces 

- Collar  grid  for  adjoining  solid  surfaces 

- Special  requirement  in  TCA  30/10  WBN  case 

• Summary 


In  this  presentation,  we  will  first  put  forth  the  background  information  for 
this  work,  including  mostly  our  past  experience  in  CFD  modeling  of  the 
TCA  W/B/N  configurations  using  face-matched  grid.  Then  we  will  describe 
the  overlapped  grid  model  for  the  TCA  0/0  WBN  configuration  and  compare 
the  CFL3D  solution  obtained  from  this  model  with  the  one  from  the  previous 
face— match  model.  Based  on  the  limited  experience  in  using  the  two  methods 
in  the  TCA  WBN  applications,  a general  comparison  of  the  two  methods  is 
presented.  Some  of  the  technical  issues  involved  with  the  overlapped  grid 
method  will  be  discussed.  Finally,  we  will  briefly  summarize  our  effort  and 
draw  some  preliminary  conclusions. 


2834 


ASE  Technologies , Inc. 


Background 

• Experience  with  CFL3D  models  using  face-match  grid  for  TCA 
WBN  configurations  with  and  without  deflected  LE/TE  flaps 

- WBN  grid  is  preferably  based  on  WB  grid 

- Grid  generation  is  time  consuming 

- Component  grids  are  not  portable 

- Grid  quality  is  limited  by  the  face-match  requirement 

- Used  wedges  to  model  deflected  TE  flaps 

• More  efficient  CFD  modeling  procedure  is  desired 

• Alternative  grid  approach:  Overlapped  grid  (MAGGIE) 


Over  the  past  two  years,  we  have  gained  a lot  of  experience  in  the  CFD 
modeling  of  HSCT  W/B/N  high  lift  configurations  using  face-matched  grid. 
Though  the  face-matched  models  have  performed  well,  the  modeling  procedure 
is  less  efficient  primarily  for  the  following  reasons.  (1)  The  W/B/N  grid  is  based 
on  the  W/B  grid  for  solution  consistency.  But  for  the  N/D  grid  to  face-match 
with  the  grid  surfaces  in  the  existing  W/B  model,  it  introduces  constraints  which 
make  the  grid  generation  process  more  time  consuming.  (2)  The  component 
grid  for  nacelles  and  diverters  can  not  be  easily  used  in  other  planforms.  (3) 
Highly  swept  wing  causes  skewness  in  the  grid.  And  interpolation  across  curved 
interfaces  introduces  numerical  errors.  (4)  Wedges  are  necessary  to  model 
the  deflected  TE  flaps  in  the  W/B/N  model,  which  cause  local  instability  in 
the  CFL3D  solution. 

The  alternative  approach  to  the  CFD  modeling  of  W/B/N  configurations  is 
the  overlapped  grid  method.  We  need  to  determine  if  the  overlapped  grid 
will  improve  the  CFD  modeling  procedure  for  TCA  W/B/N  configurations. 
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For  the  overset  grid  model  for  the  TCA  0/0  W/B/N  configuration,  we  simply 
used  the  N/D  grid  from  the  previously  completed  face-matched  model 
and  let  them  overlap  with  the  original  TCA  0/0  W/B  grid.  No  additional  grid 
generation  was  involved  for  this  overset  grid  model.  This  model  consists  of 
14  blocks  with  4.4  million  grid  points.  It  requires  180  MW  of  memory  on  Cray 
C-90  to  run  the  CFL3D  solver  for  this  model.  The  MAGGIE  preprocessor 
in  CFL3D  was  used  to  obtain  the  interpolation  coefficients  for  grid  overlapping 
between  the  N/D  grid  and  the  W/B  grid.  It  used  15  minutes  of  CPU  time 
on  Cray  C-90  for  MAGGIE  to  complete  the  interpolations. 
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TCA  0/0  WBN  Overlapped  Model 


Force  integration  on  overlapped  viscous  surfaces 

• Outer  boundary  of  N/D  grid  on  wing  lower  surface 
matches  surface  gridlines  of  the  W/B  grid 

• Overlapped  region  from  W/B  grid  on  wing  lower  surface 
is  not  included  for  lift  and  drag  force  integration 

• Black  lines  show  the  hole  definition  in  W/B  grid 
on  wing  lower  surface 


One  of  the  technical  issues  involved  in  the  overset  grid  approach  is  the 
force  integration  in  the  overlapped  region  on  model  solid  surfaces.  The 
general  pratice  is  to  go  through  additional  post-processing  steps  to  account 
for  the  right  areas  (See  C.J.  Woan’s  report  for  more  information).  In  our 
overset  grid  model  for  the  TCA  0/0  W/B/N  configuration,  since  the  N/D  grid 
is  directly  from  the  face-matched  model,  the  outer  boundary  of  the  N/D  grid 
matches  the  surface  grid  lines  of  the  W/B  grid  on  wing  lower  surface. 
Therefore,  by  excluding  the  overlapped  segment  of  the  W/B  grid  on  wing 
lower  surface  from  force  integration,  the  correct  surface  area  is  accounted  for 
in  the  overall  lift  and  drag  calculations  in  the  overset  grid  model. 
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CFL3D  Convergence  History 

(TCA  0/0  WBN,  Mach  #=0.3,  Re=8x106,  AOA=10  degrees) 


Overlapped  model: 

• No  grid  sequencing 

• 3-level  multi-grid 

• 16  hours  on  Cray  C-90 


Face-matched  model: 

• 2-level  grid  sequencing 

• 2-level  multi-grid  on  fine  level 

• 1 8 hours  on  Cray  C-90 


The  current  version  of  CFL3D/MAGGIE  allows  multi-grid  iterations  but 
does  not  support  grid  sequencing  for  the  overlapped  grid  model.  The  1200 
fine-level  iterations  completed  for  the  TCA  0/0  W/B/N  overset  grid  model 
used  16  CPU  hours  on  Cray  C-90.  With  the  face-matched  model  for  the 
TCA  0/0  W/B/N,  the  5000  coarse  level  iterations  and  500  fine  level  iterations 
used  a total  of  1 8 CPU  hours  on  Cray  C-90.  Note  that  the  original  Baldwin- 
Lomax  turbulence  model  with  Degani-Schiff  option  was  used  in  both  models 
for  the  iterations  cited  in  the  above  figure. 
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CFL3D  Convergence  History 

(TCA  0/0  WBN,  Mach  #=0.3,  Re=8x106,  AOA=10  degrees) 

1.00 
0.75 
o'  0.50 
0.25 
0.00 

0 1000  2000  3000  4000  5000  6000 

Iterations 


The  convergence  history  for  CL  and  CD  from  the  two  TCA  0/0  W/B/N 
models  is  shown  here.  For  the  last  200  iterations  in  both  cases,  the 
peak-to-valley  variations  in  CL  and  CD  are  within  2%  of  their  mean 
values. 
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TCA  0/0  WBN  CFL3D  Solutions: 

(Mach  #=0.3,  Re=8x106,  AOA=10  degrees) 


Lift  and  Drag  Coefficients 


Face-Matched 

Model 

Overlapped 

Model 

Difference 

CL 

0.3958 

0.3768 

-0.0190 

cd 

0.0627 

0.0599 

-0.0028 

For  comparison  purposes,  the  original  Baldwin-Lomax  turbulence  model 
with  Degani-Schiff  option  was  used  for  both  models.  The  changes  in  the 
aerodynamic  performance  due  to  Yeh’s  modification  to  the  turbulence  model 
were  +0.01 58  in  CL  and  +0.0007  in  CD  in  the  face-matched  WBN  CFD  model. 


For  both  the  face-matched  and  the  overlapped  grid  models,  the  mean 
values  for  CL  and  CD  from  the  last  200  iterations  are  listed  in  the  above 
table.  Apparently  there  exist  discrepancies  between  the  two  CFL3D 
solutions  in  terms  of  the  integrated  aerodynamic  coefficients.  Further 
investigation  is  necessary  to  determine  the  cause  for  the  discrepancies. 
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This  figure  shows  the  pressure  distributions  on  the  wing  upper  surface  from 
the  CFL3D  solutions  using  the  face-matched  and  the  overlapped  models. 

In  general,  the  face-matched  model  and  overlapped  model  captured  very 
similar  overall  flow  features  for  the  TCA  0/0  W/B/N  configuration. 


Similar  resemblance  in  flow  features  captured  by  the  two  models  can  be 
seen  in  the  pressure  distribution  on  the  wing  lower  surface  shown  above. 
More  detailed  contour  plots  for  pressure  coefficients  on  wing  upper  and  lower 
surfaces  may  be  necessary  to  explain  the  discrepancies  in  CL  values 
predicted  by  the  two  models. 
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Comparison  of  Face-Matched  and  Overlapped 
Grid  Interfacing  Schemes  in  TCA  WBN  Applications 


Face-Matching 

Overlapping 

Grid  Generation 

Additional  constraints  from 
surfaces  of  existing  grid 

Constrained  by  geometry  only 

Grid  Portability 

Less  flexible 

More  flexible 

Preprocessing 

Less  computer  resource 

More  computer  resource 

Grid  Interface 

2-D 

3-D 

Model  Size 

Smaller 

Larger 

Grid  Sequencing 

Yes 

No 

Multi-Grid 

Yes 

Yes 

Force  Integration 

Completed  in  CFL3D 

Additional  post-processing 
for  overlapped  solid  surfaces 

This  table  summarizes  a general  comparison  of  the  face-matched  method 
and  the  overlapped  method  based  on  our  recent  experience  with  modeling 
TCA  WBN  high  lift  configurations  using  CFL3D/RONNIE/MAGGIE.  Overall, 
overset  grid  approach  provides  more  potential  in  improving  the  CFD  modeling 
procedure  mainly  because  it  is  more  flexible  in  generating  the  component  grid 
as  well  as  in  using  the  component  grid  for  different  planforms.  As  for  the 
pre-processing  step,  RONNIE  can  usually  be  completed  in  the  debug  queue 
on  NAS  supercomputers  while  MAGGIE  has  to  be  submitted  to  the  batch  queue. 
Also,  since  CFL3D/MAGGIE  does  not  support  grid  sequencing,  it  will  be  more 
expensive  and  time  consuming  to  debug  an  overlapped  CFD  model  in  the 
development  stage. 
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Technical  Issues  with  Overlapped  Grid 

MAGGIE  Interpolation  off  Curved  Solid  Surfaces 

Curved  Solid  Surface 


Ghost  Cells 


Boundary  Cells 


Donor  Cells 


• Grid  cells  with  high  aspect  ratio  exist  near  solid  surfaces  in  large  scale  model 

• Ghost  cells  may  exist  out  of  bounds  for  even  a slightly  curved  surface 

• Several  ghost  cells  may  receive  boundary  condition  data  from  the  same 
donor  cell  in  the  overlapping  grid  through  zeroth-order  interpolation 


One  of  the  technical  issues  with  the  overset  grid  approach  is  associated 
with  interpolation  from  donor  cells  to  ghost  cells  at  the  outer  boundary  of 
the  component  grid.  In  large  scale  CFD  models  such  as  for  the  TCA  W/B/N 
configurations,  grid  cells  with  very  high  aspect  ratios  exist  near  solid  surfaces. 
When  the  solid  surfaces  are  even  slightly  curved,  the  ghost  cells  can  easily 
go  out  of  bounds  as  shown  in  the  above  figure.  In  this  case,  several  ghost 
cells  may  receive  boundary  condition  data  from  the  same  donor  cell  through 
zeroth-order  interpolation,  which  will  lose  the  gradient  information  in  the 
boundary  layer  profile  and  thus  affect  the  numerical  solution  in  the  component 
grid  as  well  as  the  overall  solution. 
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Technical  Issues  with  Overlapped  Grid 

Collar  Grid  for  Adjoining  Solid  Surfaces 


Method  1 

- Two  different  grid  faces 
for  adjoining  surfaces 

- More  flexible 


Method  2 

- One  grid  face  for 
adjoining  surfaces 

- Less  flexible 


Another  technical  issue  with  the  overset  grid  approach  is  the  usage  of  collar 
grid  for  adjoining  solid  surfaces.  In  the  literature,  the  collar  grid  described 
as  method  2 in  the  above  figure  is  often  used  and  has  been  validated. 
However,  sometimes  it  is  difficult  to  fit  this  type  of  collar  grid  in  the  actual 
configuration,  for  example,  in  the  leading  edge  region  of  nacelles  and 
diverters  adjoining  wing  lower  surface.  The  collar  grid  described  as 
method  1 in  the  above  figure  offers  more  flexibility  but  it  needs  to  be 
validated. 
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Technical  Challenge  in  Overlapped  Grid 
for  TCA  30/10  WBN  Model 


(Grid  shown  are  from  TCA  30/10  WBN  face-matched  grid) 


Finally,  the  region  between  the  nacelles/diverters  and  the  deflected  TE 
flaps  still  poses  a technical  challenge  for  overset  grid  approach  due  to 
the  close  proximity  of  the  deflected  TE  flaps  to  the  nacelles  and  diverters. 
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Summary 

• Developed  an  overlapped  grid  CFD  model  for  TCA  0/0 
Wing/Body/Nacelle  high  lift  configuration. 

• Obtained  converged  CFL3D  solution  using  the  overlapped  grid 
for  TCA  0/0  WBN  at  Mach  No.=0.3,  Re=8x106,  AOA=10  degrees. 

- Face-matched  and  overlapped  grids  yield  similar  solutions. 

- Discrepancy  exists  between  the  integrated  aerodynamic 
coefficients  using  different  grid  interfacing  schemes. 

Further  investigation  is  required. 

• Evaluation  of  the  alternative  grid  interfacing  schemes  reveals: 

- Overlapped  grid  method  is  more  adaptable  to  configuration  changes. 

- Less  time-consuming  grid  generation  with  overlapped  grid  approach. 

- Technical  challenges  still  exist  for  deflected  TE  flaps  configurations. 

- Overlapped  grid  will  be  more  efficient  for  TCA  WBN  configurations. 


In  summary,  our  recent  work  is  directed  towards  the  evaluation  of  overset  grid 
method  for  the  CFD  modeling  of  TCA  W/B/N  high  lift  configurations.  With  the  first 
overset  grid  model  for  the  TCA  0/0  W/B/N  configuration,  we  obtained  a 
converged  CFL3D  solution  which  compares  well  with  the  CFL3D  solution 
using  the  face-matched  model  in  terms  of  overall  flow  features.  Discrepancies 
exist  between  the  integrated  aerodynamic  coefficients  from  the  two  CFL3D 
solutions  and  we  will  further  investigate  the  reasons  for  the  discrepancies. 

Based  on  our  extensive  experience  with  face-matched  grid  interfacing  method 
and  the  limited  experience  with  the  overset  grid  interfacing  method,  we  conclude 
that  overset  grid  approach  offers  more  potential  in  improving  the  CFD  modeling 
procedure  for  the  TCA  W/B/N  high  lift  configurations. 
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